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Abstract 

The analysis of the experimental data of Crystal Barrel Collaboration on the 
pp annihilation in flight with the production of mesons in the final state resulted 
in a discovery of a large number of mesons over the region 1900-2400 MeV, thus 
allowing us to systematize quark-antiquark states in the (n,M 2 ) and (J,M 2 ) 
planes, where n and J are radial quantum number and spin of the meson with 
the mass M. The data point to meson trajectories in these planes being approx- 
imately linear, with a universal slope. The sector of scalar mesons is discussed 
in more detail, where, on the basis of the recent X-matrix analysis, the nonet 
classification of quark-antiquark states was performed: in the region below 2000 
MeV, two scalar nonets are fixed, that is, the basic one and the nonet of the 
first radial excitation. In the scalar sector, there are two states with the isospin 
1 = 0, which are extra ones with respect to the quark-antiquark classification, 
i.e. exotic states: the broad resonance /o(1200 — 1600) and the light cx-meson. 
The ratios of coupling constants for hadronic decays to the states irir, KK, ijr], rjrj' 
point to the gluonium nature of the broad state /o(1200 — 1600). 
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1 Introduction 



To understand the structure of strong interactions at low and moderate energies is one 
of the highest-priority problems of the modern particle physics. In the last decades, 
great efforts have been paid to develop the strong- interaction theory, the strong QCD, 
with a considerable progress though without determinate breakthrough. One may 
believe that, in part, this was due to a poor knowledge, up to recent time, of the 
experimental situation in meson sector. 

To enlighten experimental situation, a number of experiments had been undertaken 
during 1990-2000 directed purposefully to the search for new meson resonances, and a 
detailed investigation of mesons discovered before was carried out. 

The Crystal Barrel Collaboration accumulated one of the richest statistics on the 
pp and np annihilation reactions; the group of physicists from PNPI took part in the 
analysis and interpretation of the data by Crystal Barrel Collaboration (a brief review 
of this activity is given in [1]); this survey presents the results of recent investigations. 

In 1993-1994, two groups, QMWC (London) and PNPI (St. Petersburg), were analysing 
meson spectra obtained in reactions of the pp annihilation at rest. As a result, the 
scalar-isoscalar resonances / (1370) and / (1500) have been discovered [2, 3, 4, 5]; now 
these resonances are actively discussed in connection with the glueball problem. To 
establish the systematics of quark- ant iquark states, the study of pp annihilation in 
flight carried out in 1999-2001 [6] was of a particular importance, for it allowed us to 
investigate the mass region 1900-2400 MeV: a large number of resonances discovered in 
this region made it possible to fix the qq trajectories on the (n, M 2 ) and (J, M 2 ) planes 
[1, 7]. We speak about the systematics of the qq states in Section 2 of this paper. 

Scalar mesons play crucial role for understanding the strong QCD. Starting from 
1995, the PNPI group worked upon the if-matrix fit to the wave (IJ PC = 00 ++ ) based 
on the simultaneous analysis of all data available by that time. The necessity of a 
combined analysis was dictated by the existence of large interference effects "resonance- 
background" as well as the effects associated with the resonance overlapping. In a 
situation of such a type, only a combined fitting to a large number of reactions allows 
one to hope for reliable results. In Section 3, a current understanding of scalar-meson 
sector is presented based on recent results of the if- matrix fit to meson spectra [8]. 

Previous analyses carried out in 1997-1998 [9, 10, 11] were based on the experi- 
mental data as follows: 

(1) GAMS data on the S*-wave two-meson production in the reactions np — > 7r 7r°n, 
rjrjn and rjrj'n at small nucleon momenta transferred, \t\ < 0.2 (GeV/c) 2 [12, 13]; 

(2) GAMS data on the jcit S-wave production in the reaction np — > 7r°7r°n at large 
momentum transfers squared, 0.30 < \t\ < 1.0 (GeV/c) 2 [12]; 

(3) BNL data on the reaction n~p — > KKn [14]; 

(4) CERN-Miinich data on tt+tt" -> tt+tt" [15]; 

(5) Crystal Barrel data on pp(at rest, from liquid hydrogen)^ 7T 7r 7r°, 7r°7r°?7, 7i°rjrj 
[3,5,16]. 
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Now the experimental basis has much broadened, and additional samples of data 
are included into the analysis [8] of the 00 ++ wave as follows: 

(6) Crystal Barrel data on proton-antiproton annihilation in gas: pp(at rest, from 
gaseous hydrogen)-^ 7r 7r°7r°, 7r°7r°?7 [17], 

(7) Crystal Barrel data on proton-antiproton annihilation in liquid: pp(ai rest, from 
liquid hydrogen)^ tt+tt-tt , K + K~tt°, K s K s n°, K + Ksn~ [17]; 

(8) Crystal Barrel data on neutron-antiproton annihilation in liquid deuterium: np{aX 
rest, from liquid deuterium)^ 7r°7r°7r~, 7r~7r~7r + , K s K~7i°, K s K S 7i~ [17]; 

(9) E852 Collaboration data on the tttc S-wave production in the reaction n~p — > 7r°7r°n 
at the nucleon momentum transfers squared < \t\ < 1.5 (GeV/c) 2 [18]. 

In addition to Ref. [9], the reactions of the pp annihilation in gas have been also 
included into the analysis [8]. One should keep in mind that in liquid hydrogen the 
pp annihilation is going dominantly from the S'-wave state, while in gas there is a 
considerable admixture of the P-wave, thus giving an opportunity to analyse the three- 
meson Dalitz plots in more detail. New Crystal Barrel data allowed us to study the 
two-kaon channel with a more confidence as compared to what had been done before. 
This is undoubtedly important for the conclusion about the quark-gluon content of the 
scalar-isoscalar /o- m esons under investigation. 

Experimental data of the E852 Collaboration on the reaction n~p — > 7r°7r°n at 
piab = 18 GeV/c [18] together with the GAMS data on the reaction -n~p — > 7r°7r°n at 
Piab = 38 GeV/c [12] give us a solid ground for the study of the resonances / (980) and 
/ (1300), for at large momenta transferred to the nucleon, \t\ ~ 0.5 ~ 1.5 (GeV/c) 2 , 
the production of resonances is accompanied by a small background, thus allowing us 
to fix reliably their masses and widths. This is especially important for /o(1300): in 
the compilation [19] this resonance is quoted as /o(1370), with the mass in the interval 
1200 — 1500 MeV, though experimental data favour rather definitely the mass near 
1300 MeV. 

The .fT-matrix amplitude determines both the amplitude poles (masses and widths 
of resonances) and K- matrix poles (masses of bare states). The ^-matrix poles differ 
from the amplitude poles in two items: 

(i) The states corresponding to the i^-matrix poles do not contain any component 
associated with the decay processes, i.e. the transitions into real mesons. The absence 
of a cloud of real mesons allows us to refer conventionally to these states as bare ones 
[9,11,20]. 

(ii) Due to the transition bare state(l) — > real mesons — > bare state (2) the observed 
resonances are mixtures of bare states — in the first place this is related to the fo- 
mesons. So, for quark systematics, the bare states are primary objects rather than the 
resonances. 

This can be explained with the example of the behaviour of a level in the potential 
picture. Consider a potential well and the levels which correspond to stable states. 
Then gradually we switch the decay channels on, that is, replace impenetrable wall by 
potential barrier — at the beginning this leads to a broadening of levels, but masses 
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of states remain almost the same. But when the widths and resonance positions are 
such that resonances overlap, a cardinal reconstruction of the structure of levels oc- 
curs. Namely, one of resonances accumulates the widths of its resonance-neighbours 
thus making them comparatively narrow. In this way an intensive mixing of resonance 
states takes place, while the positions of masses are shifted in a value of the order 
of the resonance width, in 100-200 MeV. The iiT-matrix amplitude allows us to trace 
the transformation of stable levels (bare states) into resonances. The ^-matrix am- 
plitude works with parameters related to stable levels. Multichannel unitarization of 
the amplitude, which is inherent in the if-matrix representation, with the account for 
analyticity, allows one to trace the transformation of stable levels into real resonances. 
A characteristic feature of the i^-matrix fit is its ability to reconstruct the picture of 
stable levels as well as re-create a real picture of complex masses and partial widths. 

The latest i^-matrix fit [8] gave us rather definite information on the resonances 
/o(980), / (1300), / (1500), / (1750) and broad state / (1200 - 1600). Relying on 
the extracted partial widths for transitions to the channels nn, KK, nrj, rjr]', one can 
analyse the quark-gluonium content of these states. In this way the following properties 
of resonances are to be formulated: 

1. / (980): This resonance is dominantly the qq state, qq = nn sirup + ss cos ip, 
with a large ss component. Assuming the glueball admixture to be not greater than 
20%, Wgiuonium < 20%, the hadronic decays give us the following constraints for mixing 
angle: —95° < ip < —40°. Rather large uncertainties in the determination of mixing 
angle are due to a high sensitivity of coupling constants to a plausible small admixture 
of the gluonium. At W g i uon i um = 0, hadronic decays provide us with ip = —67° ± 10°. 

2. / (1300) (in the compilation [19] this resonance is denoted as / (1370)): This 
resonance is the descendant of the bare qq state which is close to the flavour sin- 
glet. The resonance /o(1300) is formed due to a strong mixing with the primary 
gluonium and neighbouring qq states. The quark-antiquark content of / (1300), qq = 
nn cos ip + ss sin ip determined from the transitions /o(1300) — > nn, KK, nn strongly 
depends on the admixture of the gluonium component. At W g i uon i um < 30% the 
mixing angle changes, depending on W g i uon i um and interference sign, in the inter- 
val —45° < </?[/ (1300)] < 25°; at W g i uonium = the hadronic decays provide us 
V2[/ (1300)] = -6°±10°. 

3. /o(1500): This resonance is the descendant of a bare state with large nn compo- 
nent. Like /o(1300), the resonance /o(1500) is formed by mixing with the gluonium and 
neighbouring qq states. The quark-antiquark content, qq = nn cos ip + ss sin <p , depends 
on the admixture of the gluonium: at W g i uon i um < 30% the mixing angle changes, de- 
pending on W g i uon i um , in the interval -20° < v?[/ (1300)] < 25°; at W g i uonium = one 
has v?[/ (1500)] = 11° ±8°. 

4. /o(1750): This resonance is the descendant of the bare state belonging to the 
radial-excitation nonet 2 3 P 1 qq, the wave function of which has large ss component. 
The i^-matrix analysis permits the two solutions, with different values of the ss com- 
ponent. In the first solution, the ss component dominates; in the absence of gluonium 
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component, </?[/o(1750)] = —72° ± 5°, and if the gluonium admixture does not ex- 
ceed 30%, then —110° < (p[f (1750)] < —35°. In the second solution, in the absence 
of gluonium component, <^[/ (1750)] = —18° ± 5°, and with the gluonium admixture 
< 30%, -50° < ¥>[/ (1750)] < 10°. 

5. /o(1200 — 1600): The broad state is the descendant of the primary glueball. The 
analysis of hadronic decays of this resonance confirmed its glueball nature: the glueball 
descendant has the quark- ant iquark component (qq) g iuebaii = (uu + dd+ y/Xss) /y/2 + A 
[21], which is defined by the probability to create new gg-pairs by the gluon field 
uu : dd : ss = 1 : 1 : A, the suppression parameter for the production of strange quarks 
being in the interval A ~ 0.5 — 0.8 [22, 23]. In terms of mixing angle for the nn and ss 
components, this means that the glueball descendant must have fgiuebaii — 27° — 33°, 
and just these values for <^[/o(1200 — 1600)] have been obtained in all variants of the 
i^-matrix fit [8]. Such a value of mixing angle, <y?[/o(1200 — 1600)] = 30° ± 3°, does not 
allow one to determine the admixture of the qq component in the broad state. This 
is due to the fact that the (qq) g i U ebaU and gluonium components are coupled to the 
channels nn, KK, nn and nn' in equal proportions. 

Referring to the broad state / (1200 — 1600) as a resonance needs some comments. 
The observed spectra in the channels nn, KK, rjr), t\t( demand to introduce a broad 
bump, and it occurred that this bump behaves in a universal way, thus enabling us 
to describe it as a resonance state. Indeed, a characteristic feature of the resonance is 
the factorization property: the resonance amplitude may be represented in the form 
gi n (s — M 2 Y x g out , and universal constants g in and g out depend on the sort of initial 
and final states only. The description of a large number of reactions in [8] agrees 
well with the factorization property of the broad-state amplitude. A large width of 
the /o(1200 — 1600) does not allow us to determine reliably its mass, but due to a 
strong production of the /o(1200 — 1600) in a large number of reactions we can, with a 
confidence, to find out the ratio of couplings to the channels nn, KK, 7)7), 7)7]' , that is, 
to define its content in terms of the qq and gluonium states. 

The i^-matrix [8] analysis does not point determinedly to the existence of the light 
a- meson which is actively discussed at present (e.g. see [24] and references therein), 
in particular in connection with the recently reported signals in the nn spectra of the 
decays D + -> n+n+n- (pole in the amplitude at M = (480 ± 40) - i(160 ± 30) 
MeV [25]), J/V -> nnu (pole at M = (390 1^) - «(141 ±1) MeV [26]), r -> nnnv 
(pole at M ~ 555 — «270 MeV [27]). Possible explanation of this discrepancy may 
consist in a strong suppression of the light a-meson production in the pp annihilation 
processes, like pp — > nnn, though, as one may think, there is no visible reason for 
such a suppression. Recall that the statistics in the Crystal Barrel reactions is by 
two orders of magnitude larger than in [25, 26, 27]. Alternative explanation can be 
associated with a restricted applicability of the ^-matrix approach at small invariant 
energy squared s: the i^-matrix amplitude does not describe properly the left-hand 
singularities associated with meson exchanges in the crossing channels. So one may 
suppose that the i^-matrix analysis does not reconstruct analytical amplitude correctly 
at Res < Ami, i- e - at (ReM) 2 - (ImM) 2 < 4m 2 , (Re M) 2 - (ImM) 2 < 4m 2 . In 
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numerous analyses, including those carried out in the dispersion relation technique 
where left-hand cuts can be accounted for in one way or another, the pole ascribed 
to the light (T-meson occurred just at Res < 4m^, e.g. see [28, 29, 30, 31]. In this 
connection let us emphasize that the dispersion-relation N/D analysis [31], where the 
low-energy 00 ++ amplitude was sewed with the ^-matrix amplitude [9], leads to the 
pole at M ~ 430 - i325 MeV. 

Therefore, the .fT-matrix amplitude analysis tells us that in the scalar-isoscalar 
sector there are two states, the broad resonance /o(1200 — 1600) and the light sigma- 
meson f (300 — 500), which are extra ones for the gg-systematics of mesons. In the final 
Section 4, the status of these possible exotic states is discussed as well as two more 
candidates for exotics — 7r 2 (1880) and the (J PC = 1~ + ) state (recall that gg-system 
cannot have such quantum numbers). 

2 Experimental data of the Crystal Barrel Collab- 
oration and systematics of mesons states 

In 1989-1997, the Crystal Barrel Collaboration studied the reactions of the pp-annihilation 
at rest and in flight at LEAR (CERN), that resulted in the accumulation of a huge 
statistics on multi-meson states. However, initially the analysis of experimental data 
had not been carried out in due course. The first attempt to analyse pion spectra in 
the reaction pp (at rest) — > 7r 7r 7r° based on a simplified isobar model led to a wrong 
identification of the peak near 1500 MeV; the peak had been identified as the tensor 
meson AX(1515) [32], and at the same time no scalar states were seen in the irir spectra 
over the range 1200-1700 MeV. 

In reactions such as pp (at rest) — > 7r 7r°7r° the three-particle interaction charac- 
teristics reveal themselves in full scale, that should be accounted for in the analysis 
of meson spectra. The method based on the extraction of leading singularities (pole, 
square root, logarithmic ones, and so on) in the production amplitude of a few parti- 
cles was developed in the papers [33, 34, 35]. This method had been applied with the 
purpose to re-analyse the Crystal Barrel data for the reaction pp (at rest) — > 7r 7r°7r°, 
and in the work performed together with the Collaboration members, the resonances 
/o(1370) and / (1500) had been discovered in the region 1300-1500 MeV [2, 3, 4]. The 
data on the reactions pp — > ir°7r°ri and pp — > 7i°rjrj involved into combined analysis en- 
abled us to discover a group of tensor and scalar- isovector resonances [5]. The masses 
and widths of these resonances obtained in the latest analysis are shown in Table 1, 
and the analysed reactions are quoted in Table 2. 

A combined fit to data collected in various reactions and experiments is a specific 
feature of the method used in the papers [8, 9, 10]. The matter is that, because of the 
presence of a considerable reaction background, the resonance does not always reveal 
itself as a peak in meson spectrum: due to interference effects the resonance may appear 
as a dip in the spectrum, or it may be seen as a "shoulder". Including into analysis 
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the GAMS data on the two-meson spectra in the reactions it p — > 7r°7r°n, n p — > 777771, 
7T~p -> 7777'n [12, 13], BNL data on 7r~p -> K + K~n [14] and CERN-Munich data on 
7T _ p — > 7r + 7r _ n [15] allowed us not only to determine reliably the ratios of different 
two-meson yields (that is rather important for the resonance classification), but also 
to conclude that the peak in the 7777-spectrum in the reaction -n~p — > rjrjn, which was 
previously claimed to be the resonance (2(1590) [13], appeared in fact as a result of the 
interference of the broad state / (1200 — 1600) and / (1500) resonance [36]. 

Table 1: Resonances discovered in the analysis of Crystal Barrel data. Masses and 
widths are presented in accordance with the latest combined analysis [1,8] 



Resonance 


jGjPC 


Mass, 
MeV 


Width, 
MeV 


/o(1300) (or / (1370) [21]) 

/o(1500) 

/ (1200- 1600) 

a (1450) 

a 2 (1660) 

/ 2 (1565) 


0+0++ 
0+0++ 
0+0++ 
1-0++ 
1-2++ 
0+2++ 


1310 ± 20 
1495 ± 6 
1400 ± 200 
1520 ± 25 
16701$ 
1580 ± 6 


280 ± 30 
126 ±5 
1200 ± 400 
240 ± 20 
310 ±40 
160 ±20 
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Table 2: The list of reactions used in the combined analysis [1,8] 



Crystal Barrel 
reactions 



Two-particle data 



Reaction 



Collaboration 



pp ■ 
pp ■ 
pp ■ 
pp ■ 
pn 
pn 
pp ■ 
pp ■ 
pp ■ 
pn 
pn 



ttVtt 

777T 7r 

777771"° 

7r + 7r~7r 







7T + 7T 7T 

7r~7r°7r 

ix°K s K s 

7T°K + K- 

n'KsKs 
ir°K-K s 



IT p — > 7T 7T ' U 

n~p — > 7r 7r°n 

7T _ p — > 777771 
7T~J> — > 7/77'n 

7T _ p — > KKn 
ir~p — > 7r 7r°n 



CERN- Miinich 

GAMS 

GAMS 

GAMS 

BNL 

E852 
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Table 3: Resonances discovered in the analysis of the pp-annihilation reactions in flight 
[6]. One star means that the resonance had been observed in one reaction only or it 
reveals itself poorly. Two stars mean that the resonance was seen in two reactions or 
in one where its contribution dominates. Three stars mark well-established resonances 
by using several reactions. 



Resonance 




Mass, MeV 


Width, MeV 


Status of the state 


7T 


i-o-+ 


2070 ± 35 


310 ± 80 


* 


7T 


i-o-+ 


2360 ± 30 


300 ± 80 


* 


<2l 


1-1++ 


2270 ± 50 


300 ± 70 


* 


7T2 


1-2-+ 


2005 ± 20 


210 ±40 


* 


7T2 


1-2-+ 


2245 ± 60 


320 ± 60 


* 


a 2 


1-2++ 


1950 ± 40 


180 ± 40 


** 


a 2 


1-2++ 


2030 ± 20 


205 ± 30 




^* z 


1-2++ 


2175+?° 


310 ± 60 


* 


^"Z 


1"2++ 


2255 ± 20 


230 ± 15 


*** 


«3 


1-3++ 


2030 ± 20 


150 ±20 




CI3 


1-3++ 


2275 ± 40 


150 ± 20 


* 


CI4 


1-4++ 


2005 ± 30 


180 ± 30 


*** 


CI4 


I-4++ 


2255 ± 40 


330 ± 70 




IT4 


1-4-+ 


2250 ± 15 


215 ± 25 


** 


fn 


0+0++ 


2105 ± 15 


200 ± 25 




fn 
J u 


0+0++ 


2320 ± 30 


175 ± 45 


* 


h 


0+2++ 


1920 ± 40 


260 ± 40 




fa 


0+2++ 


2020 ± 30 


275 ± 35 




fa 


0+2++ 


2240 ± 30 


245 ± 45 




fa 


0+2++ 


2300 ± 35 


290 ± 50 




fi 


O+4++ 


2020 ± 25 


170 ± 20 


*** 


fi 


O+4++ 


2300 ± 25 


280 ± 50 




no 

iz 


0+2-+ 


2030 ± 40 


190 ± 40 




iz 


0+2-+ 


2250 ± 40 


270 ± 40 


*** 


UJ 


o~i — 


2150 ± 20 


235 ± 30 




LO 


o~i — 


2295 ± 50 


380 ± 60 


* 


UJ 2 


0-2— 


1975 ± 20 


175 ± 25 


* 


U 2 


0-2— 


2195 ±30 


225 ± 40 


* 


UJ 3 


O-3— 


1960 ± 30 


165 ±30 


** 


hi 


0-1+- 


2000 ± 20 


205 ± 20 




hi 


o-i+- 


2270 ± 15 


175 ± 30 




Pi 


1+1— 


1980 ± 30 


165 ±30 




92 


1+2- 


1940 ± 40 


155 ±40 


* 


92 


1+2— 


2225 ±35 


335 ± 75 


* 


93 


1+3— 


1980 ± 15 


175 ± 20 


** 


93 


1+3— 


2260 ± 20 


200 ± 30 


* 


94 


1+4— 


2240 ± 25 


210 ±40 




bi 


1+1+- 


1970 ± 40 


215 ±60 




h 


1+1+- 


2210 ± 50 


275 ± 45 


* 


h 


1+3+- 


2020 ± 15 


110 ± 20 


** 


b 3 


1+3+- 


2245 ± 50 


350 ± 80 


* 
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The data by the Crystal Barrel Collaboration discussed above were obtained from 
the pp-annihilation at rest. In addition, the Crystal Barrel Collaboration has a huge 
statistics for the events of the pp-annihilation in flight, the antiproton momentum 
covering the range 600-1900 MeV. After the Crystal Barrel Collaboration stopped its 
activity in 1999, these data were provided to the PNPI group for further processing 
and analysis of spectra. In 1999-2001, together with the colleagues from Queen Mary 
and Westfield College (London) and Rutherford-Appleton Laboratory, the PNPI group 
analysed these data. More than thirty resonances have been discovered in the mass 
range 1900-2400 MeV [6]; they are shown in Table 3. The discovery of these resonances 
allowed us to establish systematics of meson gg-states on the (n, M 2 )- and (J, M 2 )- 
planes [1, 7] (n is radial quantum number of the gg-state with mass M and J is its 
spin) . 



2.1 Systematics of meson states 

In Figs. 1 and 2, the trajectories on the (n, M 2 )- and (7, J PC )-planes are shown for 
the states with negative and positive charge parities: 

C = - : M11+-), 6 3 (13 + -), M01+-), P(H~), p 3 (13— ), 

u/<f>(01—), w 3 (03— ) ; (1) 

C = + : tt(10- + ), 7T 2 (12-+), 7r 4 (14-+), 77(00"+), ^(02"+), 
a (10 ++ ), ai(ll ++ ), a 2 (12++), a 3 (13 ++ ), a 4 (14 ++ ), 
/o(00 ++ ), / 2 (02 ++ ). (2) 

In terms of gg-states, the mesons of the nonets n 2S+1 Lj fill in the following (n, M 2 )- 
trajectories for M < 2400 MeV: 

'So ^ tt(10- + ), ^(00- + ) ; (3) 

3 S, ^ ), w(01— )/0(Ol— ) ; 

1 P 1 -> 6i(ll + -), /ii(01 + ~) ; 

3 Pj ^ aj(U ++ ), fj(0J ++ ), J = 0,1, 2; 

l D 2 -> tt 2 (12-+), 77 2 (02-+) ; 

3 Dj -> pj(lJ— ), uj(0J—)/<l>j(pJ—), J =1,2,3; 

X F 3 ^ 6 3 (13 + -), /i 3 (03 + -) ; 

3 F 7 ^ aj(lJ++), /j(0J++), J = 2, 3, 4. 



Different orbital momenta can form the trajectories with the same J, namely, J = 
L ± 1. Therefore, the number of such trajectories doubles; these states are (J, 1 ), 
(J, 2 ++ ), and so on. Isoscalar states have two flavour components nn = (uu + dd)/\/2 
and ss, that also results in a doubling of trajectories such as ??(00~ + ), / (00 ++ ), and 
so on. 
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The trajectories with negative charge parity, C — —, can be defined practically 
unambigously (in Figs. 1 and 2, black circles stand for the observed states [6, 8, 19], 
while open circles mark the states predicted by trajectories). The trajectories are linear 
with a good accuracy: 

M 2 ~ M 2 + (n - l)/i 2 , (4) 

where M is the mass of the basic meson with n — 1, and the slope parameter is about 
(j 2 ~ 1.3 GeV 2 . The trajectory slopes for £»i(ll H ) and 6 3 (13 H ) are slightly lower: for 
them p 2 ~ 1.2 GeV 2 . 

In the sector with C = +, the states ttj belong definitely to linear trajectories with 
p 2 ~ 1.2 GeV 2 , with an exception for 7r(140) that is not a surprise, for the pion is 
rather specific state. The sector of the aj-states with J = 0, 1,2,3,4 demonstrates a 
neat set of linear trajectories, with the slopes p 2 ~ 1.15 — 1.20 GeV 2 ; the same slope is 
seen for the f 2 - and / 4 -trajectories. 

For the / - m esons the trajectory slope is p 2 ~ 1.3 GeV 2 . Let us emphasize that 
two states do not lay on linear gg-trajectories, namely, the light a- meson /o(300 — 500) 
[19] and the broad state /o(1200 — 1600), this latter has been fixed by the if-matrix 
analysis [8, 9, 10]. 

The picture of the state location on the (n, M 2 )-plot is complemented by trajec- 
tories on the (J, M 2 )-plots: they are shown in Fig. 3. To draw the daughter (J, M 2 ) 
trajectories, it was important that leading meson trajectories (tt, p, , ai, a 2 and P') were 
well-known from the analysis of hadronic diffract ive scatterings at pi a b ~ 5 — 50 GeV/ c. 

Pion trajectories, leading and daughter ones, are linear with a good accuracy, see 
Fig. 3. The other leading trajectories, p,rj, a 1: a 2 , ji or P', can be also considered as 
linear, with a good accuracy: 

a x (M 2 ) ~ a x (0) + a' x (0)M 2 . (5) 

The parameters of linear trajectories determined by masses of the gg-states are as 
follows: 

«,,(()) ~ -0.015 , a4(0) ~ 0.72 GeV" 2 ; 

a p (0) ~ 0.50 , a' p (0) ~ 0.83 GeV" 2 ; 

^(0) ~ -0.24 , a' v (0) ~ 0.80 GeV" 2 ; 

a Ol (0) ~ -0.10 , a^(0) ~ 0.72 GeV~ 2 ; 

a a2 (0)~0.45, «a 2 (0) ~ 0.91 GeV -2 ; 

a P .(0) ~ 0.71 , a' P ,(0) ~ 0.83 GeV -2 . (6) 

The slopes of the ct' x (0) trajectories are approximately the same. The inverse slope 
value, l/a' x (0) ~ 1.25 ± 0.15 GeV 2 , is of the order of the slope p 2 for the trajectories 
on the (n, M 2 )-plane: 
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The daughter trajectories for n (Fig. 3a), a\ (Fig. 3b), p (Fig. 3c), 02 (Fig. 3d) and 77 
(Fig. 3e) are determined unambigously. On the P'-trajectories, leading and daughter 
ones (Fig. 3f), there is no room for the states /o(300 — 500) and / (1200 — 1600): this 
fact stresses once again that these states are superfluous for the gg-systematics, i.e. 
they are exotic states. 

2.2 Exotics — the states which do not belong to the quark- 
antiquark trajectories 

As is seen, there are two states in the scalar-isoscalar sector which do not belong either 
to the (n, M 2 ) or (J, M 2 ) trajectories; these are the broad state / (1200 — 1600), which 
was found in the if-matrix analysis of the 00 ++ -wave [8, 9, 10], and the light cr-meson. 

The problem of the light cr-meson is discussed rather long ago [19, 28, 29, 30, 31], 
yet, until now, there is no common opinion about the existence of this state. 

Recently there appeared indications that a signal from cr-meson has been seen in 
the reactions D + — > 7r + 7r + 7r~ [25], J ftp — > it-klu [26] and r — > ixixixv [27]. Nevertheless, 
one cannot state that it was a reliable experimental identification of this state, be- 
cause in the three-particle spectra, in the region of small 7T7r-masses, an enhancement 
of the spectra can occur, due to both reflected signals from other channels and the 
rescattering effects associated with the resonance production in other channels (tri- 
angle singularity effect [34, 35]). To single out these effects rather large statistics is 
necessary, comparable with that for the reactions investigated by Crystal Barrel; how- 
ever, in the works [25, 26, 27] such a level of statistics had not been reached. Let 
us stress that the effects of reflected signals and those of triangle singularities were 
seen, when the Crysrtal Barrel reactions have been analysed, and for certain reac- 
tions they occurred to be rather important. Indeed, a correct identification of / (1500) 
in the reactions pp — > 7r 7r°7r°, 7i f]f] [2, 3, 4] became plausible after accounting for 
the interference of the decay /o(1500) — > 7T7t, 7777 and reflected resonance signals from 
other channels; the triangle-diagram effects were also studied by analysing the reac- 
tions pp — > 7r°7T 7r , 7r°?7?7, see [2, 4], but it occurred that they do not affect the spectra 
directly, and they may be effectively described by introducing complex-valued param- 
eters for resonance production. Still, in certain cases the triangle diagrams cannot 
be taken into account in such a simple way: another reaction investigated by Crys- 
tal Barrel, that is, pp — > r/7r + 7r~7r + 7r~ [37], where singular terms are located near the 
physical region, provided us with an example of strongly affected spectra. One can 
expect noticeable corrections coming from triangle diagramms also in the reactions 
D + — > 7r + 7r + 7r~, J/ip — > mvu and r — > mrnu, for the effect of triangle diagrams is 
enhanced just near the small Tin- masses, for more detail see [34, 35]. 

It is necessary to underline that in the reactions measured by Crystal Barrel the 
bumps in the 7T7r-spectra are observed, which may be considered, owing to incorrect 
treatment, as an indication to the light cr-meson. The reaction dp(at rest) — > 7r°7r°7r~ 
can serve us as an example (see [8], Fig. 5): in the 7r°7r°-spectrum, there is an obvious 
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enhancement over the phase space in the region M n o n o ~ 400 — 500 MeV. However, the 
Dalitz-plot analysis proved that this enhancement is due to the reflected signal from 
p(1450) — > 7T~7r°. The existence of a reflected signal in the D + — > 7r + 7r + 7r~ is also 
possible [38]. 

To summarize, for the experimental discovery of the light a-meson the data are 
needed, which would exceed the present statistics by one-two orders of value, i.e. it 
must be the statistics comparable with that of the Crystal Barrel reactions. 

In the paper [39], the observation of the resonance 7^(1880) was reported; this 
state does not belong to (n, M 2 )- and (J, M 2 )-trajectories. Being an extra state for 
the gg-systematics, the 7r 2 (1880) can be the hybrid: the ggg-system. 

3 Scalar-meson sector 

In this Section, the results of the i^-matrix analyses performed in [8, 9] for the scalar 
sector are subsequently presented. The sector of scalar mesons is of a particular interest; 
a variety of opinions exist concerning the properties of states which belong to this 
sector, e.g., see [11, 40, 41, 42, 43, 44, 45]. The latest analysis [8] is the most detailed 
investigation of the 00 ++ wave, where the available data have been used in a full scale. 

3.1 The if-matrix analysis of the (IJ PC = 00++)-wave 

In a set of papers [8, 9, 10, 20], the i^-matrix analysis of the waves IJ PC = 00 ++ , 
10 ++ , 02 ++ , 12 ++ had been carried out in the mass range 280-1900 MeV. For these 
states, the masses and widths of resonances had been found. In the scalar-isoscalar 
sector the following states are seen (see Fig. 4): 

00++ : / (980), / (1300), / (1500) , 

/ (1200- 1600), / (1750) . (8) 

The / (980) is a well-known resonance, its properties and its nature are intensively 
discussed during several decades. The / (1300) resonance is denoted in the compilation 
[19] as / (1370), however its mass following the most accurate determination, as is 
stressed above, is near 1300 MeV — so the notation / (1300) is used here. The /o(1500) 
resonance had been discovered in [2, 3, 4], now it is a well-established state. A few years 
ago, there existed a strong belief that in the region around 1700 MeV a comparatively 
narrow state, /j(1710), was present, with J = or 2. The ^-matrix analysis [8, 9, 
10, 20] points to the existence of / (1750), with the width r ~ 140 - 300 MeV: the 
uncertainty in the definition of the width of / (1750) is due to a bad knowledge of the 
71717171 channel in this mass range and, correspondingly, with two available solutions, 
with T ~ 140 MeV and T ~ 300 MeV. 

The broad state / (1200 - 1600), with a half-width 500 - 900 MeV, is definitely 
needed for the fT-matrix analysis. In the paper [9] this state had been denoted as 



14 



/o(1530l25o) : a large error in the definition of the mass is due to the remoteness of the 
pole from the real axis (physical region) as well as to the existence of several solutions 
given by the i^-matrix analysis. 

A large number of reactions, which were succefully described with the / (1200 — 
1600), proved the valididty of a factorization inherent in the resonance amplitude: near 
the pole the amplitude is gi n (s — M 2 )~ x g ouU where the universal coupling constants, 
gin and gout, depend on the type of the intial and final states only. A strong production 
of the /o(1200 — 1600) in various processes allows one to fix reliably these coupling 
constants. 

In the scalar-isovector sector, the i^-matrix analysis points to the presence of two 
resonances: 

10++ : a (980), a (1520) , (9) 
while in the tensor-meson sector the following states are seen: 

12 ++ : a 2 (1320), a 2 (1660) ; 

02++ : / 2 (1270), / 2 (1525), / 2 (1580) . (10) 



3.1.1 fT-matrix amplitude 

In the .fT-matrix analysis, the fitting parameters are the i^-matrix elements which are 
represented as the sums of pole terms g^gj^ j{p? n — s) (s = M 2 is the invariant energy 
squared of mesons, s = M 2 ) and a smooth s-dependent term f a b(s). Namely, 

(n)>) 

K ab = E^ + J^), (11) 

where fitting parameters are g^\ p n and f ab , and the indices a, b refer to the reaction 
channels: /o — > vr7r, KK, rjrj, rjr)' and tttttttt. The i^-matrix poles are not the amplitude 
poles, i.e. they do not correspond to real states. The amplitude in the ^-matrix 
approach is written as follows: 

t K 

T^ifx ' (12) 

where p is the diagonal matrix of phase spaces for the processes under consideration, 
p = diag (pi(s), p 2 (s), ...). Because of that, the amplitude poles correspond to the zeros 
of the determinant: 



det 



1-ipK = , (13) 



while the i^-matrix poles respond to the states with switched-off decay channels. The 
states associated with the ii'-matrix poles do not contain a cloud of real mesons which 
appear during resonance decay: this circumastance allows one to call the fT-matrix 
poles the "bare states" [9, 11, 20, 36]. 
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3.1.2 Partial amplitude for the 00 ++ wave: unitarity, analyticity and the 
problem of the left-hand cut 

The K-matrix representation of partial amplitude takes account of the well-known 
fact that at low and moderate energies the inelastic processes are dominantly two- 
particle ones. Being applied to the 00 ++ wave this means that, along with the elastic 
scattering nn — > nix (threshold at y/s = 280 MeV), we have the transitions nn — > KK 
(threshold at y/s = 2m K ), nir — > r\r\ and nn — > 7777' (thresholds at y/s = 2m v and 
y/s = m ri + m v i). Besides, in the 00 ++ -wave at y/s > 1300 — 1400 MeV, a considerable 
four-pion production is observed, but this process can be also treated, with a good 
accuracy, as a formation of two p's or two effective a-mesons: tttv — > pp — > tttttttt and 
7T7T — > a a — > 7T7r7r7r. The amplitude in the if-matrix representation takes into account 
correctly both unitarity and threshold singularities for the two-meson processes. In 
this way, the K-matrix representation of the amplitude serves us as a right frame for 
the correct reconstruction of the amplitude above the im threshold (in our case, over 
the range 280-1900 MeV). 

In the X-matrix amplitude, the threshold singularities are taken into account by 
treating the phase spaces as analytical functions of the total energy squared s: above 
the threshold, the two-particle phase space can be represented in a standard form 
p a (s) = \Jl — (mi a + m 2(l ) 2 / s, and below threshold it should be expressed through 

analytical continuation: p a (s) = i^J (m\ a + m 2a ) 2 /s — 1. The phase space of the two- 
resonance state, of the type of pp or era, may be expressed in the form which reproduces 
correctly the threshold singularities, namely, threshold singularity at y/s = 4m w and the 
singularity in the lower complex half-plane s (on the third sheet) at y/s = 1(m p — iF p /2), 
that is related to the production of pp, see [8, 9, 11] for more detail. 

The singularities, which are not explicitly taken into account within the X-matrix 
approach, are the so-called left-hand side singularities of the partial amplitude. These 
singularities are due to the exchange of particles in the crossing channels (i.e. in the t- 
and li-channels) , they determine the interaction forces. In the 7T7r-scattering amplitude, 
the nearest left singularity is located at s = 0; this singularity is associated with the 
two-pion exchanges in the t- and -u-channels. A large contribution is provided by the 
p-meson exchange, that leads to the logarithmic singularity of the partial amplitude at 
s = Ami ~ m p — — 0.5GeV 2 . In this very region, the two-pion exchange in the 00 ++ 
wave, that corresponds to the effective a-meson, contributes significantly too. The 
contribution of tensor-meson exchanges is important in the region s ~ — 1.5 GeV 2 . 

The problem of a correct account for left-hand singularities becomes even more 
complicated, for the contributions from various exchanges cancel each other at small s 
to a great extent [46], so the contribution from the left-hand cut depends strongly on 
the details of the t- and w-channel exchange mechanism, in particular on the structure 
of the p, a, ji — > 7r7r form factors. 

Because of uncertainties in choosing the interaction forces, it would be reasonable 
not to inflict any hypothesis about the left-hand cut for the 7T7r-scattering amplitude, 
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but to allow a freedom for the fitting procedure. The i^-matrix technique, which was 
used in [8, 9], allows one to do that: this technique may be easily generalized in a way 
that enables to take account of the left-hand amplitude singularities. This opportunity 
was used in [8, 9], according to what the data require; this is an important item, let us 
discuss it in more detail. 

By discussing the left-hand cut of the 00 ++ -amplitude, it is sufficient to restrain 
our consideration by one, the nearest, tttt channel only. In this case, partial amplitude 
can be written as follows: 

= ray • < 14 > 

where the pion loop diagram, which ensures the unitarity of the partial amplitude, is 

B(s) = J*?- N( ; s >^ . (15) 

J IT S S ?<0 

Ami 

At s > 4ra^, the imaginary part of the loop diagram is ImB(s) = p n7T (s)N(s) (half- 
residue in the pole s' = s), and the real part is the principal value of the integral (15). 
The real part, Re-B(s), does not contain the 7T7r-threshold singularity; this singularity 
is defined by the phase space factor in the imaginary part, p- irn {s). So we may re-write 
the amplitude A(s) in the form 

A(s) = , r , (16) 

where in the right-hand half-plane s > the i^-matrix block 

may have the pole singularities only, and in the left-hand one, s < 0, it contains a set 
of left-hand singularities. By singling out the left-hand side singularities in the explicit 
form, one can write the X-matrix block as follows: 

K(s) = 52-T— + f( s ) > ( 18 ) 
i mf — s 

where the pole positions, mf , are given by the equalities: 

[1 - Re B(s)] s=m 2 = , (19) 
and f(s) contains all left-hand singularities: 

- N ^ \- 9 2 i _ ~[°<M disc f(s') 

J[) ~ l-ReB(s) Zfmj-s ~ J tt s> - s ' [ UJ 

SL 
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The magnitude sl defines the location of the closest left-hand singularity; for the 
considered 7T7r-scattering amplitude, sl = 0. 

In the physical region, a suitable approximation for the amplitude is given by sub- 
stituting the spectral integral (20) with the sum: 



— cx 

/ 



dJ_ disc f(s') ^ y, f n ^ 

7T S' - S „ S + S n 

where s n > — s L . In the multichannel analyses [8, 9], the approximation had been used, 
when the left-hand cut in the ^-matrix terms K ab was fitted to the one-pole term as 
follows: 

/«*(*) - -4*" • ( 22 ) 
s + s 

The parameter s m different variants of the fit [8] appeared in the interval 0.5 < 
s < 1.5 GeV 2 , that points to a large contribution of the t- and w-channel f 2 , p and 
a exchanges. The use of the two-pole approximation for the multichannel if-matrix 
deteriorates the convergence of the fitting procedure, for the available experimental 
data are insufficient to fix unambigously a large number of the left-hand-cut parameters. 

More fruitful for the reconstruction of analytical structure of the amplitude near 
s ~ is the fitting to the 00 ++ wave in the region s < 1000 MeV, when one may 
restrict oneself by the one-channel approximation. In Section 3.8, we present more 
detailed narration of the study of analytical structure of the amplitude at the region 
s ~ — 4m 2 , with the restoration of the left-hand singularities obeying equations (20) 
and (21); this topic is related to the light a-meson problem. 



3.1.3 Three-meson production in the reactions of the pp and np annihilation 

The i^-matrix representation of the amplitude may be applied to the description of the 
production of resonances in the three or more particle production. The uniformity of 
the amplitude representation is important in the combined analysis of the two-particle, 
like 7T7T — > 7T7T, KK, and multiparticle, like pp — > mrn, nrjrj, nKK, reactions. The K- 
matrix approach to multiparticle processes is based on the fact that the denominator 
of the K-matrix two-particle amplitude (12), [1 — pK} -1 , describes pair interactions of 
mesons in the final state as well. 

Let us clarify this statement using as an example the amplitude of the pp annihila- 
tion from the 1 S level: pp( 1 5 l ) — > three mesons. Let the produced mesons be labelled 
by the indices 1,2,3; then the production amplitude for the resonance with the spin 
J = in the channel (1,2) is as follows: 

A 3 (s 12 ) = Kt ompt \s l2 )[l - i Pl2 K 12 (s 12 )]-\ (23) 

where the matrix factor [1 — ipuK^^s^)]^ 1 depends on the invariant energy squared 
of the mesons 1 and 2 and it coincides with matrix factor of the two-particle amplitude 
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(12). The factor omp (si 2 ) stands for the prompt production of particles 1,2 and 
resonances in this channel: 



(Kt° mpt \s 12 )) ^ = £ + Vab{si2 ) , (24) 

n 

where A^ n ) and ^> a b are the parameters of the prompt-production amplitude, and 
and \x n are defined by the two-meson scattering amplitude, see (12). 

It is also appropriate to mention that the description of the two-particle interactions 
by using the if-matrix factor (1— pK)' 1 is nothing else than a generalization of the well- 
known Watson-Migdal formula used for the two-nucleon interactions at low energies 
in nuclear reactions with multiple production of nucleons [47]. 

The whole amplitude for the production of the (J = 0)-resonances is defined by the 
sum of contributions from all channels: 



A 3 (s 12 ) + A 2 (s 13 ) + A^s) (25) 
To account for resonances with the nonzero spin J one needs to substitute in (23): 

M*u) - E4 J) («i2)4S,.^(^)4S,.^(^), (26) 
j 

where the .fT-matrix amplitude A<jP (si 2 ) is determined by the expression similar to (23), 
and Xff stands for meson states with the angular momentum J. The angular- 
momentum operators depend on the perpendicular components of meson relative mo- 
menta: and k%. Here h± 2 is the component of relative momentum of particles 1 and 
2, ki2 = (k\ — k 2 )/2, which is orthogonal to the total momentum of particles 1 and 2, 
P12 = (hi + k 2 ), namely, [k^Pii) = 0; likewise, (k^p) = 0, where p = k\ + k 2 + k 3 . 

For the lowest angular momenta, the operators Xj^ 2 [Jhj can be easily written; 
for example, for J = 1,2 we have, up to the normalization factor, X^ik^) ~ k~t 

and X ^ 2 ( k± ) ~ ( k iA ~ l k±2 9ii»2)' where is the metric tensor in the space 

orthogonal to the total momentum. The construction of operators for arbitrary J may 
be found in [48]. 

The covariant operators Xj£) are determined in the four-momentum space; 
this is a relativistic generalization of the 3-dimensional angular-momentum operators 
of Zemach [49]. The use of the 4-dimentional operators is suitable for the analysis of 
multiparticle final states, for in this case one does not need numerous Lorentz boosts, 
which are necessary in the Zemach technique or by using spherical functions. 

The amplitude expansion with respect to the states with different angular momenta 
by using relativistic covariant operators has been carried out for all the reactions pp — > 
three mesons, in the analysis of which the PNPI group took part [2, 3, 4, 5, 6, 8]. Such 
an expansion, apart from being simple and compact, as was mentioned above, has an 
advantage of taking correct account of kinematical factors, which are significant for the 
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reconstruction of the correct behaviour of the threshold singularities of multiparticle 
amplitude. In the other orbital-moment decomposition techniques this is a subject of 
a special care; in particular, in [50] a special method was suggested to take account of 
kinematical factors by using the expansion with respect to spherical functions. 

In the paper [48], the operators are also constructed for spins and total angular 
momenta in case of fermion and photon systems, which are used for the analysis of the 
reactions pp mesons and 77 — > mesons. 

The formula (23) serves us for singling out the pole singularities of the amplitude 
PP^Sq) — > threemesons, which are the leading ones. The next-to-leading logarith- 
mic singularities are related to the rescattering of mesons produced by the decaying 
resonances (that is, triangle-diagram singularities [34]). The analysis performed in 
[2, 4] showed that in the reactions pp(at rest)^ 7r°7r 7r°, ir ir n, ir°nn, to determine 
parameters of a resonance produced in the two-meson channels one may not take into 
consideration explicitly the triangle diagram singularities — it is important to account 
only for the complex-valuedness of the prompt production amplitude, that is due to 
final-state interactions (that is, the complexity of parameters A^ n ) and ip ab in (24)). 
Note that it is not a universal rule for the meson production processes in the pp an- 
nihilation — for example, in the reaction pp — > nii + 7i^7i + 7i^ [37], the explicit form of 
the triangle singularity is important. 

A complete account for unitarity and analyticity in the three-meson production 
amplitude is related to the consideration of a full set of meson interactions in the final 
state. In the reactions pp(ai rest)^ nnn, nun, nnn, nKK, the connected system of 
the dispersion-relation N/D equations, with all pair-meson equations accounted for, 
has been written in [51]. The three-meson production amplitudes being related to each 
other by the N/D equations leave less freedom for the fitting than formula (23) and, 
in principle, they provide more information about meson-meson amplitudes. However, 
the fitting on the basis of the N/ D equations is much more complicated procedure than 
the X-matrix analysis. 

3.1.4 Peripheral two-meson production in meson— nucleon collisions at high 
energies 

The two-meson production reactions up — > nun, KKn, nnn, nn'n at high energies 
and small momentum transfers to the nucleon, t, provide us with a direct informa- 
tion about the amplitudes -kit — > -kit, KK, nn, nn', for at \t\ < 0.2(GeV/c) 2 the 
reggeized n exchange dominates the production reactions. At larger \t\, the change of 
the regime occurs: at \t\ > 0.2(GeV/c) 2 a significant contribution of other reggeons 
become plausible (ai-exchange, daughter 71- and ai-exchanges). Despite a vague knowl- 
edge of the exchange structure, the study of the two-meson production processes at 
\t\ ~ 0.5 — 1.5 (GeV/c) 2 looks rather attractive, for at such momentum transfers the 
broad resonance vanishes, so the production of the / (980) and / (1300) appears practi- 
cally without background, that is important for finding parameters of these resonances. 
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The amplitude of the peripheral production of two S-w&ve mesons reads: 



(ip N R ip N )R(s 7TN ,t)K 7rR ( t ) 1-pK , (27) 

where the factor (iPnOriPn) stands for the reggeon-nucleon vertex, and Or is the spin 
operator; R(s w n, t) is the reggeon propagator depending on the total energy squared of 
colliding particles, s^n, and momentum transfer squared t, and the factor K n R^[l — 
ipK\~ x is related to the block of the two- meson production — one may find the detailed 
description of the amplitude (27) in [8, 52]. 

The factor K nR ( t )[l — ipK]" 1 describes the transitions TrR(t) — > tttt, KK, 7777, 777/: 
the block K^r^ is associated with the prompt meson production , and [1 — ipK]^ 1 
is a standard factor for meson rescattrings, see (12). The prompt-production block is 
parametrized as follows: 



where G^(t) is the /o production vertex, and f^R^ stands for the background pro- 
duction of mesons, while the parameters and p n are the same as in the transition 
amplitude txtx — > tttt, KK, rjr), rjr]', see (12). 

At the early stage of the analysis of the S'-wave two-pion production, n~p — > 
(iriT)sn, the mechanism of the reggeized-7r exchange was suggested both at small and 
moderately large \t\ [53]. The change of the regime at \t\ ~ 0.2 — 0.4 (GeV/c) 2 has been 
described by including the effective 7rP-exchange (P is the pomeron): the amplitude of 
the 7rP-exchange has another sign than the 7r-exchange, thus leading to the amplitude 
zero and, correspondingly, to a zero in the t-distribution at \t\ ~ 0.3 (GeV/c) 2 , where 
the regime changes. The filling-in of the dips in the t-distributions at \t\ ~ 0.3 (GeV/c) 2 
can appear due to exchanges of other reggeons, namely, reggeized ax-exchange and 
contributions from the daughter trajectories, ^(daughter) and ai^aughter) ■ Practically, all 
these terms do not interfere in the reaction 7i~p — > (7171)3 n: the tx- and ax-exchanges 
do not interfere due to different spin structures in reggeon-nucleon vertices, (0^ ~ (t± 
and ai ~ <T||, where cr ± and cry are the transverse and longitudinal components of the 
Pauli matrix operating in the spin space of the nucleon). At the same time the contri- 
butions of leading and daughter trajectories do not interfere practically due to a phase 
shift in reggeon propagators (for example, the propagator of the leading 7r-trajectory 
may be considered, with a good accuracy, as real magnitude, while the propagator of 
the 7T( rfaus?itef ,)-trajectory is nearly imaginary), see [52] for more detail. 

The ai exchange for the reaction np — > (7m) s n was considered in [10, 54]. In 
[10], the calculations of the mi spectra were performed both with the ai exchange and 
within an effective n exchange. As a result, it occurred that the a\ exchange affects 
weakly the parameters of the /o(980) and /o(1300): the matter is that the tttt spectra 
measured by GAMS group [12] were averaged over large ^-intervals, thence the details 
of the t-distributions are not significant for fixing resonance parameters. 
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New data on the production of the nn system at \t\ < 1.5 (GeV/c) 2 [18] gave rise 
to the discussion about the role of the t exchange in the definition of parameters of the 
/o(980) and / (1300). The data of the E852 Collaboration obtained at lower energies 
(piab = 18 GeV/c) as compared to the GAMS energy (pi ab = 38 GeV/c) point definitely 
to the fact that the description of the peripheral pion production n~p — > (7nr)sn in 
terms of the leading t exchanges, tt and a\, is not complete: the alteration of spectra 
with energy in the region \t\ ~ 0.3 — 0.4(GeV/c) 2 and M nn ~ 1300 MeV proved a 
considerable weight of daughter trajectories: it (daughter) and/or a l{daughter) . 

Combined analysis of the (irn)s spectra by GAMS [12] and E852 [18] at \t\ < 
1.5(GeV/c) 2 has been performed in [52]. The analysis [52] showed that, though the 
data do not allow us to find out the t exchange mechanism unambigously, this cir- 
cumstance affects weakly the definition of parameters of the produced /o(980) and 
/o(1300) resonances: in all variants of the fit (including various combinations of the 
^(leading), ^(leading), ^(daughter), ai(daughter) exchanges, the consideration of the effective 
irP and a\P exchanges at \t\ > 0.2 (GeV/c) 2 , with P being the pomeron, or the Orear 
mechanism [55]) the parameters of / (980) and / (1300) are almost the same. This cir- 
cumstance allows us to restrain ourselves in the K-matris fit [8] by the two variants of 
the t-channel mechanism, namely, n (leading), am eadin g), n (daughter) or ^(leading), ai(ieading), 

^l(daughter) ■ 

An opposite statement, according to which the choice of the t-exchange mecha- 
nism influences signficantly the definition of the parameters of / (980) and / (1300), 
was claimed in [56] (though without performing the fitting to resonance parameters). 
However, the analysis [52] does not confirm the statements of [56] . 



3.2 Classification of scalar bare states 

The systematics of scalar bare states has been carried out in [57], where bare .fT-mesons 
were found which are needed to fix two quark- ant iquark nonets, l 3 Poqq and 2 3 Poqq. 
The qq nonet contains two scalar-isoscalar states, /g are (l) and fl are {2), scalar-isovector 
meson a b Q Te and scalar kaon K h Q Te '. The decay couplings to pseudoscalar mesons for 
these four states, 

/o 6are (l), /o are (2) - mcKKwrrf , 
a b are -> 7TT], KK , 

K b are -> tcK, rjK , (29) 

are determined, in the leading-order terms of the 1/iV-expansion [58], by three param- 
eters only. They are the common decay coupling g, parameter A for the probability 
to produce strange quarks in the decay process (in the limit of a precise SU(3) flavour 
symmetry we have A = 1) and mixing angle for ss and nn = (uu + dd)/\/2 compo- 
nents in the / - mesons ; i> flavour (/o) — nncosip + sssirnpi. For the nonet partners, 

M are (i)] - M are m = 90°. 
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The rigid constraints on the decay couplings of bare states (29) are imposed by the 
quark-combinatorics relations. The rules of quark combinatorics were first suggested 
for the high-energy hadron production [59] and then extended for hadronic J/\I> decays 
[60]. The quark-combinatorics relations were used for the decay couplings of the scalar- 
isoscalar states in the analysis of the quark-gluonium content of resonances in [61] and 
later on in the K- matrix analyses [8, 9, 10, 11, 20, 57]. 

These constraints being imposed on the decay couplings (29) provide us with an 
opportunity to fix unambigously the states belonging to the basic nonet [57], also see 
[8,9]: 

l 3 P qq : / bare (700 ± 100), / 6are (1220 ± 40) , 

a^ re (960±30), ^ 6are (l220^ ) , (30) 

as well as the mixing angle for / o bare (700) and / 6are (1220), 

M are {700)] = -70° ± 10° , (31) 
V?[/o are (1220)] = 20° ±10°, 

To establish the nonet of the first radial excitation, 2 3 P qq, is more complicated task. 
The i^-matrix analysis gives us two scalar-isoscalar bare states at 1200-1650 MeV, 
/ 6are (1230 ± 40) and / bare (1580 ± 40), whose decay couplings (29) obey the relations 
appropriate to the glueball. The matter is that the relations between couplings for 
the glueball decay, glueball — > nn, KK, rjrj, rjr]', and the decay couplings of the quark- 
antiquark flavour singlet, (qq) singlet tttt , K K , rjr] , rjr]' , are almost the same (they are 
exactly the same in the limit A = 1). Because of that, it is impossible, by using hadronic 
decay couplings only, to distinguish between the glueball and flavour singlet. 

Systematics of bare gg-states on the (n, M 2 )-plane helps us to resolve the dilemma 
which one of these states is the glueball. This systematics (which is discussed in more 
detail in Section 3.8) definitely tells us that the state /Q are (1580 ±40) not being on the 
gg-trajectory is an extra one, so furthermore we accept that this state is the glueball. 
Indeed, 

0++ glueball : / feare (1580 ± 40) . (32) 

Lattice calculations [62] agree with this statement: gluodynamical glueball should be 
in the mass range 1550-1750 MeV. 

Having accepted the /g are (1580 ± 40) to be non-gg-state, we construct the nonet 
2 3 P qq in a unique way: 

2 3 P gg : /o are (1230±40), / 6are (1800 ± 30) , 

a^ are (1650±50), K b are (l885±fg ) . (33) 

The kaonic bare states were defined in the i^-matrix analysis of the 7ri^-spectrum 
performed in [57] on the basis of data [63]. This analysis provided us with a few 
solutions. In (30) and (33) the values of iT bare (1220l^ ) and i^^ e (1885±fg ) were 
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used, which were the average ones for these solutions: a large error in the determination 
of masses of bare states resulted from the variaty of mass values coming from different 
solutions. 

After switching on the decay channels, the bare states (30), (32), (33) turn into real 
resonances. For scalar-isoscalar states we have, after the decay onset, the transforma- 
tion as follows: 

/o(980) , 
/o(1300) , 
/o(1500) , 
/o(1200 - 1600) , 

/o(1750) . (34) 

This evolution of states is illustrated by Fig. 5, where the shift of amplitude poles into 
comlex plane is shown depending on gradual onset of the decay channels. Technically, 
it is not difficult to switch on/off the decay channels for the i^-matrix amplitude: one 
should substitute in the i^-matrix elements (11): 

9a n) - Ux)9 ( a n) , fat - tf(x)U , (35) 

where the parameter-functions for switching on/off the decay channels, £ n (^) and £/(^), 
satisfy the following constraints: £ n (0) = = an d £n(l) = — 1) an d x varies 
in the interval < x < 1. Then, at x — 0, the amplitude A turns into the K- matrix, 
A(x — > 0) — > K, and the amplitude poles occur on the real axis, that corresponds to 
the stable /Q are -states. At x = 1, we deal with real resonance; varying x from x = to 
x = 1 we observe the shift of poles into the complex M-plane. 

The x-dependence of £ n (x) and (,f(x) is governed by the dynamics of the state 
mixing, and the ^-matrix solution does not clear up this dynamics. In [64], the mixing 
has been modelled within the two-component approach for the decay processes, when 
the resonance widths are due to the transitions fo — > qq and f — > gg. This model 
proved that just the glueball accumulated the widths of neighbouring qq states. The 
dominant accumulation of widths by the glueball occurs by virtue of two reasons. First, 
mutual transitions qq-mesons <-> glueball are not suppressed within the 1/N expansion 
rules, e.g. see [11]. Second, the orthogonality of the (/(/-states suppresses direct mixing 
of the gg-mesons. 

3.3 Ovelapping /o- resonances in the region 

1200-1700 MeV: the accumulation of widths of quark- 
ant iquark states by the glueball 

The formation of the broad state is not an accidental phenomenon. The broad state ap- 
peared as a result of mixture of resonances, due to the transitions /b(l) - ► r ^ol mesons - 



/ 6are (700 ± 100) 
/ kire (1220±40) 
/ feare (1230±40) 
/ feare (1580±40) 
/ 6are (1800±40) 
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/o(2). Such transitions in case of overlapping resonances, result in a specific phe- 
nomenon: when several resonances with common decay channels overlap, one of them 
accumulates the widths of neighbouring resonances. So a broad resonance appears 
together with several narrow ones. 

Initially, the effect of accumulation of widths had been discussed in nuclear physics 
[65, 66, 67]. Concerning the scalar 00 ++ -mesons, the accumulation of widths by one of 
overlapping resonances was observed in [20, 36, 68]. In [68, 69], the following scheme 
has been suggested: the broad state /o(1200 — 1600) is the glueball descendant; this 
state was formed because of the glueball mixing with neighbouring gg-states, that 
was accompanied by the accumulation of widths of neighbouring states by the glueball 
descendant. As a result, comparatively narrow states / (1300), / (1500), / (1750) have 
had considerable admixtures of the glueball component, while the broad state got a 
large gg-component. The quark- ant iquark component in the glueball should be close 
to the flavour singlet, namely [21]: 

(qq) glueball = (uU + dd + V\ss) / \Jl + A, (36) 

where the parameter A is nearly the same as that in the decay processes, A ~ 0.5 — 0.8. 

In meson physics, the accumulation of widths can play a decisive role for the destiny 
of exotic states which are beyond the gg systematics. Indeed, the exotic state, after 
appearing in a set of gg-states, creates a group of overlapping resonances. If the 
transition of the exotic state into gg-state is not suppressed (as in case of the glueball, 
where, according to the 1/iV-expansion rules [58], the transition glueball — > qq-meson is 
allowed in the leading-order terms [11, 70]), then just the exotic meson has an advantage 
to accumulate the widths: the wave functions of neighbouring gg-states are orthogonal 
to each other but not to the exotic state. Therefore, the existence of the broad state 
together with comparatively narrow ones must serve as a signature of exotics in this 
mass region [71]. 

The broad state may play rather constructive role in the formation of the confine- 
ment barrier. The broad state, after accumulating the widths of resonances-neighbours 
in the mass scale, plays the role of a locking state. The evaluation of radii of the broad 
state / (1200 — 1600) and two narrow neighbouring ones, / (980) and / (1300), which 
was performed in [10, 72], tells us that the radius of the broad state is considerably 
larger than the radii of / (980) and / (1300): this fact agrees well with the assumption 
that /o(1200 — 1600) plays the role of a locking state for its resonance-neighbours. Re- 
cent measurement of the t-distributions in the reaction n~n — > (mv) s n [18] confirmed 
relatively large magnitude of the broad state radius: with the increase of \t\ (momen- 
tum square transferred to the resonance), the broad state vanishes much faster than 
/o(980) and / (1300). 
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Table 4: Coupling constants squared (in GeV 2 ) of scalar-isoscalar resonances decaying 
to the hadronic channels 7T7r, KK, rjr], rjr)' and mnm for different K-matrix solutions. 



Pole position 


-nit KK r/r] rjr]' nimn 


Solution 


/o(980) 
1031 - i32 
1020 - i35 


0.056 0.130 0.067 0.004 
0.054 0.117 0.139 0.004 


I 

II 


/o(1300) 
1306 - il47 
1325 - il70 


0.036 0.009 0.006 0.004 0.093 
0.053 0.003 0.007 0.013 0.226 


I 

II 


/o(1500) 

1489 - i51 

1490 - i60 


0.014 0.006 0.003 0.001 0.038 
0.018 0.007 0.003 0.003 0.076 


I 

II 


/o(1750) 
1732 - i72 
1740 - ilQO 


0.013 0.062 0.002 0.032 0.002 
0.089 0.002 0.009 0.035 0.168 


I 

II 


/o(1200 - 1600) 
1480-ilOOO 
1450 - i800 


0.364 0.265 0.150 0.052 0.524 
0.179 0.204 0.046 0.005 0.686 


I 

II 



3.4 Evolution of couplings of the 00 ++ -states to channels tttt, 

ttttttjt, KK, TjTj r]r]' with the onset of decay processes 

The if-matrix analysis does not allow one to determine partial widths of the Jo- 
resonances directly. To find out partial widths for the decays f — > im, ixixixix, KK, 
rjr), rjr)' it is necessary to calculate the residues of the amplitude poles corresponding to 
resonances. Near the resonance, the transition amplitude a — > b (indices a and b stand 
for the resonance channels nn, KK, rjr], r\r( ', numi) takes the form: 

(n) ("■) 

A ab c e^+^) + Bab . (37) 

The pole position defines the mass and width of the resonance /z n = M n — i(T n /2), 
and the real- valued coupling constants to channels, and g^, allow us to find 
partial widths of resonances. Factorized residues are complex-valued magnitudes, their 
complexity is determined in (37) by the phases 6^ and 9^; in (37), a smooth, non- 
pole, term B ab is also written. 

On the basis of the results of the K- matrix fit [8], the decay coupling constants 
0<W were calculated for / (980), / (1300), / (1500), / (1750) and the broad state 
/o(1200 — 1600), and the comparison was done for the obtained values with corre- 
sponding couplings of bare states, which are predecessors of the resonances under 
discussion. 

The values of couplings calculated in [8] are shown in Table 4. The comparison 
reveals a significant difference between the decay couplings for bare states and their 
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descendant -resonances. This undoubtedly proves a strong effect of the mixing of qq- 
states with the glueball: the real resonance is a mixture of these states. 

Figure 6 demonstrates the evolution of coupling constants at the onset of the decay 
channels: following [73], relative changes of coupling constants are shown for / (980), 
/o(1300), /o(1500) and / (1750) after switching on/off the decay channels (recall that 
the value x = corresponds to the amplitude poles on the real axis and the value x — 1 
stands for the resonance observed experimentally). 

Let us bring our attention to a rapid increase of the coupling constant /o — > KK 
at the trajectory /Q are (700) — /o(980) in the region x ~ 0.8 — 1.0, where 7 2 (x = 
1.0) — 7 2 (x = 0.8) ~ 0.2, see Fig. 6a. Actually this increase is the upper limit of 
the possible admixure of the long-range KK component in the / (980): it cannot be 
greater than 20%. 

3.5 Evaluation of the glueball component in the resonances 



/o(980), / (1300), / (1500), / (1750) and broad state / (1200- 



1600) based on the analysis of hadronic-decay channels 

The evolution of couplings observed in the transition from bare states to real resonances 
is due to the mixture of the gg-states with the glueball that is a consequence of the 
transitions /o(l) — > real mesons — > /o(2). One can evaluate the quark- ant iquark and 



glueball components in / (980), / (1300), / (1500), / (1750) and / (1200- 1600) using 



the rules of quark combinatorics for the decay couplings. 

In the leading-order terms of the 1/iV-expansion the vertices of hadronic decays of 
the resonances are defined by planar diagrams, the examples of the planar diagram 
for the decay of gg-state and the glueball into two mesons are presented in Fig. 7a,b, 
respectively. In the course of the gg-state decay, the gluons produce a new gg-pair; in 
the glueball decay, a subsequent production of two pairs occurs. 

For the decay couplings squared for / — > tut, KK, rjr), 777/, the quark-combinatorics 
rules, in case when the /o state is the mixture of the quarkonium and gluonium com- 
ponents, give us [8, 73]: 





2 



COS(f + sin ip sin 0) + 



+ A 



(cos 2 + A sin 2 0) 




(38) 
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The terms proportional to g stand for the transitions qq — > two mesons, while those 
with G respond to transitions glueball — > two mesons. Accordingly, g 2 and G 2 are 
proportional to the probability to find in the considered / - meson the quark-antiquark 
and glueball componets. Recall that the angle cp stands for the content of the qq- 
component in the decaying state, qq = cos ipnfi + sin Lpss, and the angle G for the 
contents of n and n' mesons: n = cos Qnn — sin ss and n' = sin 6 nn + cos 9 ss; we 
use = 38° [74]. 

One may believe that the decay of the glueball is going in two steps: initially, one 
qq pair is produced, then with the production of the next qq pair a fusion of quarks 
into mesons occurs. Therefore, at the intermediate stage of the f decay, we deal with 
a mean quantity of the quark-antiquark component, (qq), which later on turns into 
hadrons. The equation (38), under the condition G — 0, defines the content of this 
intermediate state (qq) = nncos((p) + ss sin (<£>). 

As was said above, the i^-matrix analysis [8] gave us two Solutions, I and II, 
which differ mainly by the parameters of the resonance /o(1750). Fitting to the decay 
couplings squared for these Solutions leads to the values of (</?) as follows: 

Solution I: 



/o(980) 
/o(1300) 

Broad state / (1200 - 1600) 
/o(1500) 
/o(1750) 



(v?)~-68°, A ~ 0.5 -1.0, 

(ip) ~ (-3°) - 4° , A -0.5 -0.9, 

(ip) ~ 27° , A ~ 0.54 , 

(ip) ~ 12° - 19° , A ~ 0.5 -1.0, 

(v?)~-72°, A ~ 0.5 -0.7, 



(39) 



Solution II: 



/o(980) : (<p) ~ -67° , A ~ 0.6 - 1.0 , (40) 

/o(1300): (if) ~ (-16°) - (-13°) , A ~ 0.5 - 0.6 , 

Broad state / (1200 - 1600) : (<p) ~ 33° , A ~ 0.85 , 

/o(1500) : (v?)~2 -ll°, A ~ 0.6 -1.0, 

/o(1750): ((^)~-18 , A-0.5. 



In both Solutions, the average values of mixing angle for / (980) coincide with one 
another, with a good accuracy. Still, one should inderline that equations (38) allow us to 
get one more magnitude for the mixing angle of the /o(980), namely, (</?[/o(980)]) ~ 40°. 
The fact that the decay constants for / (980) — > nn and / (980) — > KK accept the 
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solution with (</?[/ (980)]) ~ 40° was underlined in [75]. However, this value of mixing 
angle does not suit the classification of bare states provided by the if-matrix solutions. 
Indeed, the / (980) is the descendant of the bare state /o are (700±100) which is close to 
the flavour octet. The evolution of coupling constants (see Fig. 6) tells us that / (980) 
by its content remains close to its predecessor, /Q are (700 ± 100). Because of that, in 
Eqs. (39) and (40) only solutions with (y>[/ (980)]) ~ -67.5° are kept. 

The values of average mixing angles for / (1300) are stable negative for both Solu- 
tions I and II, they differ slightly, so we may accept (</?[/ (1300)]) = —10° ± 6°. 

Also the mean mixing angle for the /o(1500) does not differ noticeably for Solutions 
I and II, so we may adopt (y>[/ (1500)]) = 11° ± 8°. 

For the /o(1750), Solutions I and II provide different mean values of mixing an- 
gle. In Solution I, the resonance /o(1750) is dominantly ss system; correspondingly, 
(<£>[/o(1750)]) = —72° ± 5°. In Solution II, the absolute value of mixing angle is much 
less, (v?[/o(1750)]> = -18°±5°. 

For the broad state, both Solutions give proximate values of mixing angle, namely, 
(y?[/o(1200 - 1600)]) = 30° ± 3°. This magnitude favours the opinion that the broad 
state can be treated as the glueball descendant, because such a value of the mean 
mixing angle corresponds to (p g i U ebaii = sin" 1 \J\/(2 + A) at A ~ 0.50 — 0.85. 

Let us emphasize that the coupling magnitudes for the /cr r esonances found in [8] do 
not provide us any alternative variants for the glueball descendant. Indeed, the value 
which is the closest to the (p S ingiet is the limit value of the mean angle for / (1500) 
in Solution I: (</?[/ (1500)]) = 19°. Such a magnitude being used for the definition of 
^glueball corresponds to A = 0.24, but this suppression parameter is much lower than 
those observed in other processes: for the decaying processes we have A = 0.6 ± 0.2 
[11, 22], while for the high-energy multiparticle production it is A ~ 0.5 [23]. In this 
way, the quark combinatorics points to the one candidate only, that is, the broad state 
/ (1200 — 1600); we come back to this important statement later on. 

Generally, the formulae (38) allow us to find ip as a function of the coupling constant 
ratio G/g for the decays glueball — > mesons and qq-state — > mesons. The results of 
the fit for / (980), / (1300), /o(1500), / (1750) and the broad state / (1200 - 1600) 
are shown in Fig. 8. 

First, consider the results for / (980), / (1300), / (1500), / (1750) shown in Fig. 
8a for Solution I and in Fig. 8c for Solution II. The bunches of curves in the (<p, G/g)- 
plane demonstrate correlations between mixing angle values and the G/g ratios, for 
which the description of couplings given in Table 4 is satisfactory. A vague dissipation 
of curves, in particular noticeable for / (1300) and / (1500), is due to the uncertainty 
of A in Eqs. (39) and (40). 

The correlation curves in Fig. 8a, c allow one to see, on a qualitative level, to 
what extent the admixture of the gluonium component in / (980), / (1300), / (1500), 
/o(1750) affects the quark-gluonium content, qq = nn cos p> + sssin<^ determined from 
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hadronic decays. The magnitudes g 2 and G 2 are proportional to the probability to find 
out, respectively, the quarkonium and gluonium components, W qg and W g i uonium in a 
considered resonance: 



9 



qq" qq J 



G 2 = G 



gluonium^ gluonium 



(41) 



qq 



The results of the K-matrix fit obtained in [8] tell us that the coupling constants g\ 
and G 2 luonium are of the same order of magnitude (also see the discussion in [11, 64, 70]), 
therefore we accept as a qualitative estimate: 



G 2 /g 2 



W, 



gluonium / ^^qq 



(42) 



The figures 8a, c show the following permissible scale of values (p for the resonances 
/o(980), /o(1300), /o(1500), /o(1750), after mixing with the gluonium component. 

Solution I: 



gluonium [/o(980)] < 15% 
T gluonium [/o(1300)]<30% 

T gluonium 

[/o(1500)] < 30% 
T gluonium [/o(1750)] < 30% 



-93° < </?[/ (980)] < -42°, 
-25° < ^[/ (1300)] < 25° , 
-2°<^[/ (1500)]<25\ 
-112° < y?[/ (1750)] < -32° 



(43) 



Solution II: 

W gluoniu m [/ (980)] < 15% 

W^I«on<«m[/o(1300)] < 30% 

W gluoniu m [/ (1500)] < 30% 

T gluonium 

[/o(1750)] < 30% 



-90° < ¥»[/ (980)] < -43°, (44) 
-42° < y?[/ (1300)] < 10° , 
-18°<v9[/o(1500)]<23°, 
-46° < y»[/o(1750)] < 7° . 



The ^-dependence of G/g- is linear for / (980), /o(1300), / (1500), / (1750). An- 
other type of the correlation takes place for the state which is the glueball descendant: 
the correlations curves for this case form in the (</?, G/g)-plane a typical cross. Just 
this cross appeared for the broad state / (1200 — 1600) for both Solutions I and II, see 
Fig. 8b,d. 

The appearance of glueball cross by correlation curves in the (tp, G / g)-plane is due 
to the formation mechanism of the quark-antiquark component in the gluonium state: 
in the transition gg — > (qq) giuebaii the state (qq) g i U ebaii is fixed by the value of A. So the 
glueball descendant is the quarkonium-gluonium composition as follows: 

gg cos 7 + {qq) g i U ebaii sin 7 , 

where 

[qq] glueball = Tin COS (p glueball + SS Sm (p glueball , 
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and (pgiuebaii = tan 1 y A/2 ~ 26° — 33° for A ~ 0.50 — 0.85. The ratios of couplings for 
the transitions gg — > tttt, KK, 7)7), 7)7)' are the same as for the quarkonium (qq) g i ue baii ~ * 
7T7T, KK, rjr], rji]', so the study of hadronic decays only do not permit to fix the mixing 
angle 7. This property - similarity of hadronic decays for the states gg and (qq) g i U ebaii 
- implies a specific form of the correlation curve in the (if, g/G)-plane: the gluonium 
cross. Vertical component of the gluonium cross means that the glueball descendant has 
a considerable admixture of the quark- ant iquark component (qq) g i U ebaii- Horizontal line 
of the cross corresponds to the dominant gg component. The value of A which affects 
the cross-like correlation on the (if, g/G)-plane is denoted from now on as \ g i ue bau- For 
Solution I, we have X g i U ebaii = 0.55, while for Solution II X g i U ebaii = 0.85. 

At not a large shift of A from its mean value \ g i ue bau, the coupling constants 
/o(1200 — 1600) — > 7T7T, KK, r]r), r]r( can be also described, with a reasonable accu- 
racy, by Eq. (38); in this case correlation curves on the (if, <//G)-plane take the form 
of hyperbola. Shifting the value of A in |A — \ g \ U ebaii\ ~ 0.2 breaks the description of 
couplings of the broad state by formulae (38). 

The cross-type correlation on the (if, g/G^-plane in the description of coupling con- 
stants /o — > 7T7T, KK, r)T), 7)7)' by formula (38) is a characteristic signature of the glueball 
or glueball descendant. And vice versa: the absence of the cross-correlation should 
point to the quark- ant iquark nature of resonance. Therefore, the .fT-matrix analy- 
sis proves definitely that /o(1200 — 1600) is the gluonium descendant, while /o(980), 
/o(1300), /o(1500), /o(1750) cannot pretend to be the glueballs. 

The analysis proves that /o(1300), /o(1500) are dominantly the nn-systems. Still, 
in Solution II the qq component of the resonance /o(1300) may contain rather large ss 
component in the presence of the 30% gluonium admixture in this resonance. As to the 
/o(1500), the mixing angle (</?[/ (1500)]) in the qq component may reach 25° at G/g ~ 
—0.6 (Solution I) that is rather close to if g i U ebaii- However, in this case the description 
of coupling constants g\ (Table 4) is attained as an effect of the strong destructive 
interference of the amplitudes (qq) — > two pseudoscalars and gg — > two pseudoscalar s . 
This fact tells us that one cannot be tempted to interprete / (1500) as the glueball 
descendant. 



3.6 The light cr-meson: Is there a pole of the 00 ++ -wave am- 
plitude? 

Effective a-meson is needed in nuclear physics as well as in effective theories of the 
low-energy strong interactions — and such an object exists in a sense that there exists 
rather strong interaction, which is realized by the scattering phase passing through 
the value 5g = 90° at M wn ~ 600 - 1000 MeV. In the naive Breit-Wigner-resonance 
interpretation, this would correspond to an amplitude pole; but the low-energy nn 
amplitude is a result of the interplay of singular contributions of different kinds (left- 
hand cuts as well as poles located highly, /o(1400— 1600) included) , so a straightforward 
interpretation of the a-meson as a pole may fail. 
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The question is whether the u-meson exists as a pole of the 00 ++ -wave amplitude. 
(See also [76], where this problem was particularly underlined.) However, until now 
there is no definite answer to this question, though this point is crucial for meson 
systematics. 

The consideration of the partial S'-wave nn amplitude, by accounting for left singu- 
larities associated with the t- and w-channel interactions, favours the idea of the pole at 
Res ~ 4to^. The arguments are based on the analytical continuation of the K-matrix 
solution to the region s ~ — Afi 2 [31]. 

In [31], the nn- amplitude of the 00 ++ partial wave was considered in the region 
y/s < 950 MeV. The account for the left-hand singularities had been done within the 
method described in Section 3.10.2, namely, the i^-matrix amplitude was fitted in the 
form 

where /„, /, s n ,g 2 , M 2 are parameters. The left-hand cut was fitted to six pole terms; 
the pole at s = M 2 (M ~ 900 MeV) corresponds to 5$ = 90°. The fitting was performed 
to the low-energy scattering phases, 5®, at y/s < 450 MeV, and the scattering length, 
Oq. In addition at 450 < y/s < 950 MeV the value 5^ was sewn with those found in 
the X-matrix analysis [9]: from this point of view the solution found in [31] may be 
treated as analytical continuation of the K- matrix amplitude to the region s ~ — 4to^.. 
The analytical continuation of the X-matrix amplitude of such a type accompanied by 
simultaneous reconstruction of the left-hand cut contribution provided us with the 
characteristics of the amplitude as follows. The amplitude has a pole at 



the scattering length, 
and the Adler zero at 



~ 430 - i325 MeV , (46) 
a° ~ 0.22 m; 1 , (47) 



y/s ~ 50 MeV . (48) 

The errors in the definition of the pole in solution (46) are large, and unfortunately they 
are poorly controlled, for they are governed, in the main, by uncertainties when left- 
hand singularities are restored. As to experimental data, the position of pole is rather 
sensitive to the scattering length value, which in the fit [31] was taken in accordance to 
the paper [77]: a[j = 0.26 ±0.06 to" 1 . As one can see, the solution [31] requires a small 
scattering length value: Oq ~ 0.22 to" 1 . New and much more precise measurements of 
the i^ e 4-decay [78] provided a° = (0.228 ± 0.015)to^ 1 , that agrees completely with the 
value (47) obtained in [31]. Such a coincidence favours undoubtedly the pole position 
(46). 

So, the AV-D-analysis of the low-energy 7T7r-amplitude sewn with the i^-matrix 
one [9], provides us with the arguments for the existence of the light a-meson. In a 
set of papers, by modelling the left-hand cut of the 7T7r-amplitude (namely, by using 
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interaction forces or the t- and -u-channel exchanges), the light cx-meson had been also 
obtained [28, 29, 30], but the mass values are widely scattered, e.g. in [79] the pole has 
been obtained at essentially larger masses, y/s ~ 600 — 900 MeV. 

The observation of the light a-meson is aggravated by the existence of the Adler zero 
in the 7T7r-amplitude near the 7T7r-threshold. Therefore, as was emphasized above, for 
the reliable determination of the cr-meson it is necessary to study the 7T7r-production 
in the annihilation and decay reactions (of the type of pp — > unit and D — > -nun), 
where the Adler zero is absent. However, one should emphasize again, in the pp- 
annihilation, where the statistics is rather high, the a-meson is not seen, while in the 
reaction D + — > 7r + 7r~7r + [25] the statistics is not sufficient for the reliable analysis of 
the low-energy 7T7r-spectrum. 

By discussing the search for the light cr-meson in the three-particle reactions, it is 
necessary to accentuate a requirement, had it been fulfilled, one can confidently confirm 
the discovery of this resonance. I mean the rescattering of pions in the final state. The 
matter is that the fitting to the near-threshold state, that is actually the a-meson, 
should be carried out under the correct account for the unitarity at small masses of 
the 7T7r-system, while the contribution of crossing channels (for example, resonance 
production) leads to the violation of unitarity. The account for rescatterings in the 
7T7r-channel reconstructs the two-particle unitarity. Besides, the rescatterings restore 
logarithmic singularities of the three-particle amplitude, that may affect significantly 
the region of small tttt masses, in particular, in the presence of heavy resonances. The 
effects of pion rescatterings were studied in [4], where the combined analysis of the 
Dalitz-plot has been done for the reactions pp(ai rest)^ 7r°7r 7r°, 7i°r]ri — it appeared 
that they act rather strongly upon the region of small 7T7r-masses, though they are not 
important for the analysis of resonances at 1300-1500 MeV, that was the main goal of 
the investigation performed in [4]. An opposite point of view, namely, the statement 
about the possibility to single out the signal from #-meson without taking account of 
rescatterings, was claimed in [80]. 

Concerning the problem of search for the light a-meson in the three-particle reac- 
tions, an important point should be emphasized, which being fulfilled would allow us 
to speak with confidence about reliable determination of this state. This is rescattering 
of pions in the final state. The matter is that fitting to a near-threshold state, such as 
cr-meson, should be carried out with correct account for unitarity at small 7T7r-masses; 
at the same time the contribution of crossing channels (for example, the production of 
resonances) violates the unitarity. The rescatterings in the 7T7r-channel being accounted 
for restore the two-particle unitarity near the 7T7r-threshold. Besides, rescatterings lead 
to logarithmic singularities of the three-particle amplitude, which in the presence of 
heavy resonances affect significantly the low-mass region. The effect of pion rescat- 
terings were considered in [4] at simultaneous analysis of Dalitz-plots in the reactions 
pp(atrest) — > 7r 7r°7r°, ir°i]r]: it occurred that they affect strongly the low-mass region 
though being not important for the analysis of resonances at 1300-1500 MeV — the 
main scope of investigation in [4]. An opposite viewpoint, namely, one can reliably 
detect the a-meson signal without accounting for 7T7r-rescatterings, is presented in [80]. 
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3.7 Systematics of scalar states on the (n, M 2 )- and ( J, M 2 )- 
planes and the problem of basic multiplet l 3 P qq 

The figures 1,2 and 3 demonstrate linear behaviour of the gg-meson trajectories in the 
(n, M 2 )- and (J, M 2 )-planes for a variaty of states. In this way, it would be instructive 
to juxtapose the nonet classification given by the K-matrix analysis [8, 9, 57] with the 
systematics of the gg-states in the (n, M 2 )- and (J, M 2 )-planes. 

3.7.1 The fT-matrix classification of scalars and gg-trajectories in the (n, M 2 )- 
planes 

Consider the variant which follows directly from the i^-matrix calssification of bare 
states: let us accept that the light cr-meson does not reveal itself as the amplitude pole 
or, if it exists, it is an extra state for the gg-systematics. For this case, the location of 
scalar states (fo, ao, K ) on the (n, M 2 ) trajectories is shown in Fig. 9. 

Figure 9a demonstrates the trajectories for the resonance states 00++, 9 ++ and ±0+, 
while in Fig. 9b the trajectories for corresponding bare states are shown (recall that 
the doubling of the / -trajectories occurs due to the existence of two components, ss 
and nn). It is seen that the trajectory slopes for bare and real states are nearly the 
same (in Fig. 9 the trajectory slope is fx 2 = 1.3 GeV 2 ). 

The state /Q are (1580 ± 50) certainly does not belong to any of /Q are -trajectories 
shown in Fig. 9b, the trajectories of the real /o-states does not need the state /o(1200 — 
1600) too. This is natural, provided the / bare (1580 ± 50) is the glueball and / (1200 - 
1600) the glueball descendant. 

Gluodynamical lattice calculations tell us that the lightest scalar glueball is located 
in the mass region 1550-1750 [62]. There are also other arguments in favour of just 
in this mass region one may encounter the gluonium: the estimate of the mass of soft 
(or effective) gluon points to the gluon mass m g i uon ~ 700 — 1000 MeV. Experimental 
estimations of the effective-gluon mass are based on the study of hadronic spectra in 
radiative decays J/ip — > 7 + hadrons and T — > 7 + hadrons [81, 82]; they give for 
the gluon mass the value 800-1000 MeV. The close values were obtained within model 
calculations of the quark-gluon interactions in the soft region giving the value 700- 
800 MeV for the effective-gluon mass [83, 84]. The lattice calculations of the effective 
gluon agree reasonably with these estimates: accordingly, m g i uon ~ 700 MeV [85] . It 
is natural to believe that the mass of the lightest scalar glueball is of the order of the 
double effective gluon mass, M g i uebaU ~ 2m g i uon . 

One should pay attention to the fact that the resonances /o(980), /o(1500) and 
/o(1300), /o(1750) lay neatly on linear trajectories, with a slope which nearly coincides 
with slopes of isovector trajectories: pj, aj and This circumstance argues in favour 
of the admixture of the glueball components not leading to strong violation of the 
linear behaviour of trajectories, at least for the / - meson sector. 



34 



It would be pertinent to recall that in the i^-matrix analyses [8, 9] three Solutions 
have been found denoted as I, II-l and II-2. Solutions II-l and II-2 occurred to be 
close to each other in resonance characteristics; yet, in Solution II-l the /g are (1580) 
is the gg-state, and the glueball mass is around 1250 MeV. We reject this Solution, 
for it contradicts both the linearity of /o are -trajectories and gluodynamical calculation 
results. 

Therefore, from the point of view of the classification of states on the (n, M 2 )- 
trajectories, the scheme where the /Q are (700 ± 100) and its descendant, /o(980), are 
the lightest states in the l^P^qq nonet looks self-consistent. This scheme agrees with 
the coupling constants for the transitions /g are , a^"" 6 , K h Q Te — > rnr, K K , r/r] , r/rj' [8, 11, 
57]. Moreover, model calculations of masses of the gg-mesons also point to a possible 
location of the l 3 Pogg- n o n et in the range 900-1300 MeV [86]. Still, by considering the 
nonet classification, the question arises about large ss-component in the lightest qq 
state: 67% in / bare (700 ± 100) and (45 - 75)% in / (980). But it would be appropriate 
to recall that we have already faced similar situation in the pseudoscalar sector: the 
rj- meson also contain a large ss-component, about 40%. Both states, 77-meson and 
/Q are (700 ± 100), are close to the flavour octet, and strong interaction in the flavour 
octet states may cause the large ss component in light mesons. In other words, the 
value of the ss component in lightest mesons depends on the structure of the short- 
range forces, and to determine the structure of forces is the goal of the analysis and 
the motivation to perform the systematics of meson states. 

So the root of the matter is whether the a-meson, if it exists, is the standard 
gg-state, a partner to / (980). Then the nonet of scalar mesons looks as follows: 
} (300-500), / (980), a (980) and the low-lying k meson, k(800-1000). The argument 
against this scheme is the systematics of the i^-mesons in the (J, M 2 ) plane, that is 
discussed below. 

3.7.2 Systematics of kaons in the (J, M 2 ) plane 

The existence of the scalar /t-meson, with the mass 800-1000 MeV, is under discussion 
as long as the discussion of cr-meson, e.g. see [24] and references therein. The low- 
energy data on the 7rfT-scattering are poor, and this is the main reason of uncertainty 
in the /t-meson problem. The K-matrix analysis of the 7ri^-amplitude [57] does not 
supply us with unambigous answer, for there exist solutions without amplitude pole at 
800-1000 MeV but there are also solutions with the pole at ~ 1000 MeV [57]. The K- 
matrix analysis of the high-statistics Crystal Barrel data on the reactions pp — > KKn 
and np — > KKn does not require the k- meson [8] . 

Let us turn to the systematics of kaon states on the (J, M 2 ) plane: such systematics 
of kaons brings additional argumets to the discussion of whether a and k, mesons are 
the standard qq states, members of the l 3 Pogg multiplets, or not. 

In the kaonic sector, the experimental data are poor, and a complete picture of the 
kaon disposition on the trajectories with J p = 0~,2~,4~,... and J p = l + ,3 + ,5 + ,... 
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cannot be unambigously reconstructed: in Fig. 10a,d a conventional picture is pre- 
sented for the location of these trajectories. More definite is the situation in the 
sectors with J p = + ,2 + ,4 + ,... and J p = 1~,3~,5~..., see Fig. 10b, c: here we know 
the states laying on the leading and first daughter trajectories — and just on these 
trajectories the k meson should lay, provided it is the standard gg-state. 

In the sectors with J p = 0", 2", 4", . . . and J p = 1+, 3+, 5 + , . . . (Fig. 10a,d) the 
doubling of trajectories takes place because of the presence of two spin states: for 
example, two states with J p = 1 + are defined by the quantum numbers (L — 1, S — 1) 
and (L — 1, S — 0); likewise, two states J p = 2~ are formed by two sets of quantum 
numbers, (L — 2, S — 1) and (L = 2, S = 0). Figures 10a, d represent a conjectural 
disposition of trajectories and demonstrate how many unknown states (open circles) 
are present in the groups J p = 0~, 2", 4~ and J p = 1 + , 3 + , 5 + . 

More reliably determined are the trajectories for the groups with J p = + ,2 + ,4 + 
(Fig. 10b) and J p = 1~,3~,5~ (Fig. 10c), where, as was said above, the leading and 
first daughter trajectories are fixed firmly. The leading and first daughter trajectories in 
these groups are degenerate: the states fT 2 +(1430), if 4 +(2045) and K* (890), #3- (1780), 
i^ 5 -(2380) lay on the common leading trajectory ai eading (J) ~ 0.25 + 0.90 M 2 (the 
slope is in GeV units), and the states K +(U30), If 2 + (1980) and if* (1680) lay on the 
common trajectory adaughter(J) — —2.0 + 0.90 M 2 , where J p = + , 1~, 2 + , 3 _ . 

Figure 10b demonstrates also the location of the k meson: one can see that this 
state does not lay either on leading or daughter trajectories. In order to put the k 
meson onto the qq trajectory, the existence of the 2 + kaon with the mass ~ 1600 MeV 
is needed, but the analysis [63] does not point to the presence of tensor resonance 
in this mass region. The fact that there is no room for the k meson on the (J, M 2 ) 
gg-trajectories is the argument against the construction of the lowest scalar nonet by 
using the states / (300 - 500), / (980), a (980) and k(800 - 1000), that is actually an 
indication that the k meson does not exist. 

3.8 Exotic scalar states, / (1200 - 1600) and / (300 - 500) 

Two scalar states, the broad resonance / (1200 — 1600) and the light cr-meson / (300 — 
500), remain beyond the considered here qq classification that is based on the if-matrix 
analysis [8, 9]; they should be treated as exotic states. 

The K- matrix analysis gives us an unambigous interpretation of the / (1200 — 
1600). This state is the descendant of a pure glueball which, by accumulating of 
widths of neighbouring scalar-isoscalar states, turned into the broad resonance. All 
comparatively narrow resonances from the region ~ 1500 MeV (they are /o(1300), 
/o(1500) and / (1750)) fit well the gg-trajectories in the (n, M 2 ) plane, so just the 
broad resonance / (1200 — 1600) is an extra one from the point of view of the qq- 
systematics. 

The fact that the glueball turned into the broad state allows us to say that the 
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glueball is "melted". The idea that the bound state of gluons can disappear due to 
strong interaction in the soft region was brought out rather long ago [87]. However, the 
"melting" observed in the if-matrix analysis is a specific one: the transformation of 
the glueball into the broad state occurs as a result of the decay processes (transitions 
of resonances into real mesons) that takes place at large distances, of the order or 
greater than R CO nfinement- Besides, the amplitude pole associated with /o(1200 — 1600) 
did not go far away from the physical region: the resonance half-width, that defines 
the imaginary value of mass related to the pole, is of the order of 500-1000 MeV, so 
the pole occurs inside the hemicircle, where the analysis of experimental data allowed 
us to reconstruct the analytical amplitude (supposing the threshold singularities be 
correctly taken into account). 

The transition of the lightest scalar glueball into the broad resonance / (1200— 1600) 
gives rise to a number of questions. The glueball tranformation into the broad reso- 
nance, is it a unique event (in a sense, an occasional event) or is this a common 
phenomenon for exotic states? The resonance / (1200 — 1600), after having accumu- 
lated the widths of neighbouring gg-states plays, with respect to them, the role of a 
locking state, does this lead to the increase of the proper size of / (1200 — 1600)? These 
questions were put in [10, 71]; still, to have a veracious response one needs much more 
accurate data. Qualitative evaluation of the radius of /o(1200 — 1600) had been carried 
out in [10, 71, 72], on the basis of the GAMS data for the M^-distributions at different 
momenta transferred to the nucleon in the reaction n~p — > n7i n at pi a b = 38GeV/c 
[12]. According to this estimate, the broad state is much more loosely-bound system 
than its neighbours-resonances. Recent measurements performed by the E852 Collabo- 
ration [18] support the fact that comparatively narrow resonances / (980) and / 1300) 
are more compact than the broad state / (1200 — 1600); the discussion of these data 
can be found in [8, 52]. 

The situation with the light a-meson, /o(300 — 500), is less definite as concern its 
experimental status and the understanding of its nature. The nature of the a-meson, 
if it exists as an amplitude pole, is rather enigmatic. The light a-meson is hardly 
the glueball-like formation. Also it is difficult to imagine the light a-meson to be the 
standard gg-system: the arguments against such an interpretation were formulated in 
Section 3.8. 

It was suggested in [88] that the existence of the light a-meson may be due to a sin- 
gular behaviour of the forces between quark and antiquark at large distances (in quark 
models they are conventionally called "confinement forces"). The scalar confinement 
potential, which defines the spectrum of the gg-states in the region 1000-2000 MeV, 
behaves at large hadronic distances as V^ n j inement {r) ~ ar, where a ~ 0.19GeV 2 . In 
the momentum representation, such a growth of the potential is associated with the 
singular behaviour at small q: 

V^ finement (q) ~ j A • (49) 
In the colour space, the main contribution comes from the component c = 8, i.e. 



37 



the confinement forces should be the octet ones. The question that is crucial for 
the structure of cr-meson is as follows: is there a component with the colour singlet, 
VcolfinementiQ) > i n the singular potential (49)? If the singular component with c = 1 
exists, then it must reveal itself in hadronic channels as well, that is, in the 7T7r-channel. 
In hadronic channels, this singularity should not be exactly the same as in the colour 
octet ones, because the hadronic unitarization of the amplitude (which is absent in 
the channel with c = 8) should modify somehow the low-energy amplitude. One may 
believe that, as a result of the unitarization in the channel c = 1, i.e. due to the 
account for hadronic rescatterings, the singularity of V^lf inement (q) may appear in the 
7T7r-amplitude on the second sheet and one may believe that this singularity is what we 
call "the light a- meson". 

Therefore, the main question consists in the following: the V^lf inernent (q 2 ) , does 

it have the same singular behaviour as v}^j inement (q 2 )7 The observed linearity of the 
(n, M 2 )-trajectories, up to the large-mass region, M ~ 2000 — 2500 MeV, see Figs. 1 
and 2, favours the idea of the universality in the behaviour of potentials V^lj inement 

and V^j inement at large r, or small q. To see that (for example, in the process 7* — > 
qq, Fig. 11a) let us discuss the colour neutralization mechanism of outgoing quarks 
as a break of the gluonic string by newly born gg-pairs. At large distances, (2.0- 
2.5 fm), that corresponds to the formation of states with masses 2000-2500 MeV, 
two-three new gg-pairs should be formed. This follows from the value (n c h(W)) in the 
e + e~-annihilation: at W ~ 2 GeV we have (n c h) — 3. It is natural to suggest that 
a convolution of the quark-gluon combs governs the interaction forces of quarks at 
large distances, see Fig. llbc. This means that the potential V^lj inement working at 
such large distances contains two or three t-channel gg-pairs. The mechanism of the 
formation of new gg-pairs to neutralize colour charges does not have a selected colour 
component, see e.g. [89]. In this case all colour components 3 (g) 3 = 1 + 8 behave 
similarly, that is, at small q 2 the singlet and octet components of the potential are 
uniformly singular, V c ( ^ finement (q 2 ) ~ V c ( ^ finement (q 2 ) ~ 1/g 4 . This is seen in Fig. 11a: 
the quark-gluon ladder ensures the t-channel flow of colour charge C = 3, so quark- 
antiquark interaction block being the convolution of ladder diagrams 3 <8> 3 = 1 + 8 
contains two equivalent, singlet and octet, components. This points to a similarity of 
V (1) f - t and V (8) f - 

confinement confinement 

So, from the point of view of mechanisms of colour neutralization related to the ideas 
of the t-channel production of new quark-ant iquark pairs, the existence of the light 
cr-meson is well-grounded. Therefore, it is of utmost importance to discover whether 
there exist singularities at low energies in the 00 ++ -amplitude, which we call the light 
cr-meson, and what is the type of these singularities, if any. 
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4 Conclusion 



Recent analysis of meson spectra in the pp-annihilation in flight [6] resulted in the 
opportunity to systematize gg-mesons, practically leaving too small a room for spec- 
ulations about the existence of exotic states. First of all, this is related to the 00 ++ 
sector, where all comparatively narrow resonances /o(980), /b(1300), /o(1500), /o(1750) 
lay on linear gg-trajectories on the (n, M 2 )- and (J, M 2 ) planes, that does not provide 
us with a ground for hypotheses on their exotic origin. 

In quark models of 60-70's the idea of dominant LS-splitting was pushed forward 
[90, 91] - within such a structure of forces the P-wave gg-multiplet with J = is 
lighter than multiplets with J = 1,2. Recently the models became popular, where 
/o(980) and a (980) are considered as exotic states (see [92, 93, 94] and references 
therein): in this case the l 3 Pogg-nonet should be placed in the region 1300-1700 MeV. 
It is seen that the i^-matrix analysis [8,11,12,13] renders us back to a former picture: 
the l 3 Pogg multiplet is the lightest one among the P-wave gg-states. 

Beyond the 00 ++ sector, the candidate for the exotics is the 7r 2 (1880) [39] - for 
sure this state does not belong definitely to the gg-trajectory, and it could be a hybrid 
state. Of course, the existence of this resonance should be confirmed by its observation 
in other reactions. Besides, if 7T2(1880) was the qqg-staXe, there should exist the other 
nonet members near 1900 MeV in addition to those that lay on standard (n, M 2 ) 
trajectories. However, no 2 h mesons have been observed yet, which could be extra 
ones with respect to the gg classification. 

In the literature, there are indications to a plausible existence of the non-quark- 
antiquark mesons with J PC = 1~ + (see [19] and references therein), though experimen- 
tal evidence for these states is rather poor. 

In the 00 ++ meson sector, two states are beyond the gg systematics, namely, the 
light a-meson and the broad state / (1200 — 1600). 

Concerning the problem of the light a-meson, though important and having a long 
discussion history, there is no direct evidence of its existence — on the level of fixing a 
pole in the 00 ++ amplitude on the basis of the experimental data with large statistics. 
Such a situation existed a few years ago, see the discussion in [76], yet there is no 
noticeable experimental progress till now. 

Nevertheless there are reasons which make us to concern about the problem of 
the light a-meson: this is singular behaviour of the quark-ant iquark interaction block 
at small momentum transfers, 1/g 4 , in the coordinate space this corresponds to the 
linear growth of the potential, ~ ar, at large r (in quark models such a potential is 
conventionally called "confinement potential"). Linear behaviour of meson trajecto- 
ries in the (n, M 2 )-plane at large masses (Figs. 1,2) tells us that such a behaviour 
occurs actually at rather large r, up to r ~ 2.0 — 2.5 fm. Assuming that both colour 
components of the "confinement potential", octet and singlet ones, behave at large 
r similarly, V^ et (r) ~ v}irl g i e t( r ) > we S a i n singular behaviour at small masses in the 



39 



white (hadronic) channel. Hadron unitarization of the amplitude (the account for the 
7T7r-rescatterings) , which is necessary in the white channel but is absent in colour one, 
is capable to "hide" the white-channel singularity under the nn branching cut thus 
re-creating the picture of the light cr-meson. The mechanism of colour neutralization, 
where the t-channel formation of new qq pairs is of principal importance, favours the 
hypothesis about similar growth of potentials V^ ] nfinement {r) ~ v}^ finement (r) . 

As concern the broad state /o(1200 — 1600), we know neither its mass nor the width 
reliably. Still, we do know that this broad state is strongly produced in a number of 
reactions. Also we know that, in terms of the qq and gluonium components, this 
state is for sure the mixture of the gluonium with the quark- ant iquark state which 
is close, by its content, to the flavour singlet (qq) g i U ebaii = (nny/2 + ssV\)/y/2 + A 
with A ~ 0.5 — 0.8: this fact is proved rather confidently by the relations between the 
couplings / (1200 — 1600) — > nir, KK, rjr), 777/, which are reliably found from the data. 

Another reliably determined feature of the broad state /o(1200 — 1600) is that it 
is a more loosely-bound system than the surrounding / - resonances ; this fact points 
to a possibility to create it as a result of the accumulation of widths of neighbouring 
resonances. If so, and the ^-matrix analyses of the 00 ++ waves [8, 9] performed on an 
extended data basis favours this scenario, then this circumstance explains us where the 
lightest scalar glueball vanishes, though it promised to be around ~ 1500 MeV but it is 
definitely missing among the observed narrow /o- m esons. In this way, let us emphasize 
that the possibility for a new comparatively narrow 00 ++ resonance near 1000-1800 
MeV is strictly excluded by the available experimental data. 

A scenario for / (1200 — 1600) being the glueball descendant, which lost a part of its 
gluonium component due to the mixing with neighbouring resonances, was suggested 
in [68, 69], on the basis of data which were much poorer than the present ones. The 
nowaday data which are richer not only in statistics but also in a number of investigated 
reactions, fit to this picture rather well, without noticeable contradictions to it. One 
more check of the picture with / (1200 — 1600) as a glueball descendant is expected 
in the nearest future for the high-statistics radiative decays J/^f — > 7 + hadrons: in 
hadronic spectra, the production of comparatively narrow resonances /o(1300) and 
/o(1500) in the 00 ++ wave should be accompanied by a strong "background" due to 
the production of the broad state / (1200 — 1600). 

At the time being the mixing of /o- m esons at 1300-1700 MeV is intensively dis- 
cussed, see recent papers [95, 96] as an example. However, it should be specially 
stressed that the mixing considered in the X-matrix technique is in principle different 
from mixing discussed in [95, 96]. In the i^-matrix approach, the mixing occurs due 
to transitions of bare states to real mesons — in other words, due to imaginary parts 
of hadron loop diagrams. As a result of such a mixing the amplitude poles "move" in 
the complex mass plane, a loss or accumulation of widths take place, that is, mixed 
states "push" each other in the complex plane by sinking deeply into complex plane 
or springing out to real axis; as a result, real mesons- resonances are created. In this 
way the shift of real part of resonance mass may be considerable — of the order of 



40 



resonance width. The mixing discussed in the papers [95, 96] happens without decay 
processes, that is, without imaginary parts of hadronic loop diagrams, so correspond- 
ing amplitude poles are on real axis: by mixing the levels repulse each other, with the 
increase or decrease of the mass but not width; after such a mixing the bare states 
should be created but not real resonances, as was suggested in papers [95, 96]. Once 
again: in the .fT-matrix technique the formation of the observed resonances takes place 
by the onset of the decay processes only. 

In accordance with gluodynamical lattice calculations, the second scalar glueball 
should exist at 2100-2200 MeV [97]. In [98], it was suggested that / (2105) is either 
this second glueball or strongly mixed with it. However there is not enough data to 
judge about the glueball being in this mass region: the discovered resonances lay on the 
qq trajectories quite comfortably. It is rather probable that the second scalar glueball 
turned into the broad state too by mixing with the neighbouring qq mesons — in this 
case to fix it experimentally the measurement of a large variety of spectra is needed 
as well as the evaluation of ratios for the nn, KK, rjr], rjr]', rj'rj' yields in the region 
2100-2200 MeV. 
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Figure 1: Trajectories on the (n, M 2 ) planes for the states with (C = — ). Open circles stand for the 
predicted states. 
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gure 2: Trajectories on the (n,M 2 ) planes for the states with (C = 
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Figure 3: (J, M 2 ) planes for leading and daughter trajectories: a) 7r-trajectories, b) ai-trajectories, 
c) p-trajectories, d) a 2 -trajectories, e) ^-trajectories, f) P'-trajectories. 
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Figure 4: Complex M plane in the (IJ PC = 00 ++ ) sector. Dashed line encircle the part 
of the plane where the i^-matrix analysis [8] reconstructs the analytic K-matrix am- 
plitude: in this area the poles corresponding to resonances / (980), / (1300), / (1500), 
/o(1750) and the broad state /o(1200 — 1600) are located. On the border of this area 
the light a- meson denoted as /o(450) is shown (the position of pole corresponds to 
that found in the N/D method [31]). Beyond the i^-matrix analysis area, there are 
resonances / (2030), /o(2100), /o(2340) [6]. 
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Figure 6: The evolution of normalized coupling constants 7 a = g a / \jY.b 9b a ^ the onset 
of the decay channels for / (980), / (1300), / (1500), / (1750). Curves demonstrate 
the description of constants by formula (38). 
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Figure 8: Correlation curves on the (<p,G/g) and (ip,g/G) plots for the description of 
the decay couplings of resonances (Table 4) in terms of quark-combinatorics relations 
(38). a,c) Correlation curves for the gg-originated resonances: the curves with appro- 
priate A's cover strips on the (ip,G/g) plane. b,d) Correlation curves for the glueball 
descendant: the curves at appropriate A's form a cross on the (<p,g/G) plane with the 
center near ip ~ 30°, g/G ~ 0. 
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Figure 9: Linear trajectorieson the (n, M 2 ) plane for scalar resonances (a) and bare scalar states 
(b). Open circles correspond to the predicted states. 
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Figure 10: The (J P ,M 2 ) planes for kaonic sector (open circles stand for the predicted states). 
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a b 

Figure 11: a) Quark-gluonic comb produced by breaking a string by quarks flowing 
out in the process e + e~ — > 7* — > qq — > mesons, b) Convolution of the quark-gluonic 
combs — an example of diagrams describing interaction forces in the gg systems at 
r ~ 2.0 fm. 
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